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ABSTRACT 


A method to characterize particles in rocket exhaust plumes is developed. The 
particle velocity, size, and material composition are determined from crater characteristics 
resulting from impacts into aluminum and copper targets passed through the plume. The 
targets are mounted on a steel arm approximately 21 inches (53 cm) long which is rotated 
through the plume at sufficient velocity to prevent material failure resulting from thermal 
effects. A Scanning Electron Microscope (SEM) with secondary x-ray detectors is used to 
determine the particle material, and a standard optical measurement microscope is used to 
determine the crater diameter and depth. The crater diameter and depth are used, in turn, 
as inputs to a ballistics computer code to estimate the velocity and size of the particle. 
The target has a safe residence time in the plume of approximately 50 ms before reaching 
an unacceptably high temperature. The arm must reach a velocity of 104 ft/s (32 m/s) 
before entering the plume to produce the design residence time of 20 ms. The arm is 
actuated by a torsion spring with a 5-inch (13 cm) outer diameter, 0.625-inch (16 mm) 
wire diameter, and 11 coils. A prototype of the entire rocket exhaust particle impact 
characterization system (PICS) was constructed and statically tested. 
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CHAPTER 1 


INTRODUCTION 

Statement of Problem 

Presently, particles from rocket exhausts are typically collected by placing sticky 
tape far downstream of the nozzle exit to capture low velocity particles that have survived 
the severe aerothermal environment of the plume. This procedure suffers from the 
problem that particles from the ambient surroundings are entrailed in the far plume and 
captured on the tape. It is difficult to determine which particles come from the 
surroundings and which ones come from die rocket engine. This method also provides no 
information on particles that may have been consumed in the extremely high temperatures 
of the exhaust near the exit. 

Another method, used in some vertical test facilities, involves shooting a dart with 
sticky tape through the exhaust (Sambamurthi, 1996). The particles impact the tape while 
the tape is in the plume and become attached to it. Unfortunately, this method also 
collects particles from the ambient surroundings and possibly also from the ground when 
the dart lands. 

Typical rocket test facilities have the capability to determine the elemental 
composition of the exhaust gas by spectroscopy (Olive, 1988). The major deficiency of 
this method, as applied to the current problem, is its inability to determine the molecular 
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composition of particles in the exhaust. Spectroscopy can determine which basic 
elements are present but cannot determine how they are combined. It is very important to 
know which molecules are exiting the rocket since numerous molecules can contain any 
given element. 

The current problem is to determine the size, number density, and material 
composition of particles from the nozzle exhaust flow. Several significant measurement 
problems associated with the high gas temperatures and supersonic velocity of the no zzl e 
flow are identified as follows: 

1. Temperatures can be in excess of 3,000°R (1,667 K) in the nozzle exhaust flow. 
Table 1 lists melting temperatures of potential particulates. Since several of these 
materials have melting temperatures well below 3,000°R (1,667 K), it is necessary to 
capture the particles as soon as possible, or they may melt or ablate away. It is not 
possible to place the probe in the current design into the nozzle upstream of the 
nozzle exit plane because of geometry constraints and high temperatures. Therefore, 
for best results, the particles should be characterized at the exhaust plane of the 
nozzle. 

2. The gas and particle velocity ranges also present a problem. Typical exhaust velocities 
for hybrid rocket engines are approximately 8,000 ft/sec (2,438 m/sec). A particle 
number density of approximately 100 particles/ft 3 (3.5x1 O' 5 particles/cm 3 ) is the 
initial design condition. It is desired to have distinct craters on the target; therefore, 
tne aesirea number of craters is of the order of twenty. This design crater number 
density results in target residence times in the plume on the order of 20 ms. Since the 
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probe moves only through approximately 180°, this residence time presents a 
significant problem in accelerating and decelerating the probe without breaking it. 
Particle velocities in the range of 1,000 to 8,000 ft/sec (305 to 2,438 m/sec) are 
expected. Impacting particles must have enough velocity to create deep craters. It is 
possible that particles may coalesce and form large particles without sufficient 
velocity to form craters in the target. 


Table 1. Approximate Melting Temperatures of Some Potential Exhaust Particles 


Potential Particle 

Hydroxl-T erminated Polybutadiene (HTPB) 
A1 2 0 3 
Steel 

Aluminum 

Copper 

Titanium 


Approximate Melting Temperature 



800 

4,180 

3,226 

1,680 

2,442 

O CAA 


2. The high temperatures may also cause chemical reactions between particles in the 
flow and the exhaust gases. If the particles do not impact the target before these 
reactions are complete, the sample collected will not be typical of the particles in the 

exhaust flow. 

3 . The system should be inexpensive, easy to operate, and reliable. 


Background 


Bowman (1976) presents a good overview of the uses for probes in combustion 
research. He indicates that probes have been used in many research applications other 
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than sampling particles exiting a rocket engine. They are commonly used for temperature 
measurement as well as gas, two-phase, and velocity sampling. Pneumatic probes, hot- 
wire anemometers, and microphone probes may be used for measuring velocities in 
various types of flows. Probes used in combustion research are usually water-cooled. 
Cemansky (1976) has used probes to sample NO and N0 2 levels in gas turbine 
combustors with limited success. His research indicates some of the problems, such as 
quenching effects, which are inherent in gas sampling probes. Since the current work is 
primarily concerned with particle characterization techniques, these problems will not be 
discussed here. Bilger (1975) and Bennet (1976) provide good overviews of that topic. 

Since this study is concerned with the characterization of particles from rocket 
exhausts, the literature review will focus on this subject. All published work located 
during the literature review involves capture of AI 2 O 3 particles, since these particles are a 
primary constituent of solid rocket fuel. As early as 1980, the Jet Propulsion Laboratory 
at the California Institute of Technology was interested in collecting particles from the 
Space Shuttle’s solid rocket motor exhaust (Strand et al., 1980). They used quartz-crystal 
mass monitors, Nucleopore filters, sticky-tape impactors, photoelectric cloud detectors, 
and electrical mobility analyzers to measure particulate concentrations and mass 
distributions. Strand found that the sticky-tape impactors and the electrical mobility 
analyzers provided the most useful data. Cofer et al. (1985) also investigated particles in 
the exhaust of the Space Shuttle. Their experiment involved flying an airplane, equipped 
with a Teflon filter, downward in a tight spiral through the exhaust of mission STS-61A 
immediately after launch. Any particles they collected were then examined with an SEM 
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to determine their size. Particles captured were between 0.1 pm and 10 pm in diameter. 
Particles smaller than 0.1 pm were not visible to the electron microscope, and it is 
believed that particles larger than 10 pm fell to Earth before they could be collected. 

Laredo et al. (1991) investigated plume particle size measurements in solid 
propellant micromotors. Their research involved collecting particles in the chamber, 
nozzle, and plume of a subscale rocket. Particles from the chamber and nozzle were 
collected using special scattering techniques. Supersonic particles were collected by 
placing an impact probe consisting of a stainless steel wedge coated with copper into the 
flow for 0.5 seconds. Particles as large as 50 pm were common; however, some as large 
as 85 pm were found. It was found that particle breakup was dominant in the throat and 
the chamber and that collision coalescence dominated downstream of the throat. 

Kreutle (1978) determined that particles collected after solid rocket firings were 
roughly spherical. This research involved collecting particles from the surfaces of the 
rocket after a test firing. Infrared optical properties have also been employed to analyze 
particle properties (Konopka et al., 1983). Konopka gives an excellent overview of 
research to this end. However, his work resulted primarily in determining certain 
properties of AI 2 O 3 particles. 

A more recent report relating to the use of impact probes to capture particles in 
rocket e xhau sts was published by Sambamurthi (1995). His method was to launch darts 
covered with sticky tape into the far downstream region (-150 ft or 45.7 m from the 
nozzle exit plane) of a large-scale solid rocket motor. The darts were 30 to 48 inches 
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(76.2 to 121.9 cm) long and were launched from the ground with a launch pressure of 
800 psig (5.5 MPa) and at an angle of 80° with the horizontal. Many of the darts did not 
land nose forward and were unusable. The other darts were “mostly uncontaminated.” 
Sambamurthi’s research indicates that A1 2 0 3 particles in the far plume region of a full- 
scale solid rocket exhaust are approximately 8 to 1 1 pm in diameter. This report details 
an attempt to take Sambamurthi’s idea of passing a probe quickly thr ou gh a rocket 
exhaust plume to capture particles. Implementation of this idea and operating parameters 
of the detailed procedure differ significantly from Sambamurthi’s method. 

Proposed Approach 

This study addresses the problems of particle characterization by passing a probe 
through the rocket exhaust plume to encounter the particles. The high-velocity particles 
impact the relatively soft target surface of the probe and leave craters on its surface. The 
craters are then examined under an SEM equipped with secondary x-ray detectors to 
determine particle composition, and a standard optical measurement microscope is then 
used to measure crater depth and diameter. The crater depth and diameter data are used 
as inputs for a ballistics impact code to predict the velocity and size of the~ particles! 
Since this method collects only high-velocity particles, which normally do not exist in the 
ambient environment, the possibility of collecting particles from the ambient environment 
is eliminated. 

Hie probe consists of an steel arm that supports a target of al umin um and copper 
at the end of the arm. The arm and target move transverse to the flow as shown in 
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Figure 1, where the probe arm is rotating perpendicular to the plane of the paper. The 
probe is actuated with a torsion spring and released by a pneumatic clip controlled by a 
delay generator receiving its input from the rocket firing circuit. The entire apparatus is 
called the Particle Impact Characterization System (PICS). This study outlmes the design 
procedure for the PICS. The procedure includes a model of the thermodynamic 
properties of the plume, a model of the kinetic response of the particles m the plume, and 
a study of the thermal response of both the arm and target to the severe aerothermal 
environment. Also included are the kinetics calculations for the arm and the torsion 
spring design, as well as instructions on how to manufacture the PICS and the targets. 

The proposed method obtains information on particles in the high velocity, high 
temperature plume environment. The experimenter can obtain information on particle 
characteristics as a function of axial location in the plume by operating several probes at 


various axial locations. 





CHAPTER 2 


OVERVIEW OF PARTICLE IMPACT CHARACTERIZATION SYSTEM 

Probe and Probe Actuation Mechanism 

Figure 1 presents the basic configuration of the PICS. The arm is actuated by a 
torsion spring and stopped by an automobile tire, as seen in Figures 2 and 3. The PICS 
consists of the following: 

1 . the probe arm which holds the target while it moves through the plume 

2. a torsion spring which accelerates the arm 

3. a mandrel which provides support to the torsion spring when it is deflected 

4. a probe arm collar which connects the arm to a shaft and allows the arm to rotate 
around the shaft 

5. a shaft which keeps the probe arm collar and spring mandrel in proper alignment 

6. a baseplate which can be bolted to a surface for support 

7. three vertical supports which hold the shaft in place and are welded to tne oasepiate 

8. a slider which holds the arm in the cocked position (Figures 3b and 4) 

9. a pneumatic clip which releases the slider when a current is applied to an electric 

valve on the clip (Figure 4) 

10. a delay generator which activates the electric valve 

1 1 . a chain which connects the slider to the arm. 
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Figure 3. PICS - Initial and Terminal Positions 




^ ^ ^ ^ ^ ft 


Plunger 


Figure 4. Cross Section View of the Pneumatic Clip with the 
Locked in Place 


The pneumatic clip is pressurized to 70 psig (482.6 kPa) with the slider in place inside the 
clip. A chain is wrapped around the probe and spring arm. A standard trailer winch 
(designed for recreational boat trailers) may be used to cock the spring if desired. Once 
the spring has been deflected enough to produce the desired launch velocity, the cha i n is 
connected to the slider by means of an industrial-strength spring clip. At this point, the 
arm is ready to fire. When the electric valve is opened by the delay generator, pressure 
drops inside the chamber and the pin is allowed to move perpendicular to the slider. Force 
from the torsion spring pulls the slider towards the spring, pushing the pin up in the 
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process and releasing the arm. At full deflection, the spring moves through a 135° 
included angle with the horizontal and the arm moves because of its inertia during the last 
45° of travel. It is during this period that the bottom of the arm impacts the tire. The tire 
absorbs the energy of the arm and simultaneously keeps the arm from impacting anything 
rigid which could damage either the arm or the target. Inherent to this design is the fact 
that the unloaded spring can act as a stop which prevents the arm from rebounding into 
the plume as it recoils from the tire. The automobile tire was empirically determined to be 
capable of bringing the probe to rest safely. 

Target Configuration 

The basic design of the target is illustrated in Figure 5. It consists of two half- 
cylinders which form a stepped cylindrical shaft 1 inch (2.54 cm) in length. The front 
0.5 inch (1.27 cm) of the shaft is 0.5 inch (1.27 cm) in diameter, and the rear half is 
0.25 inch (0.635 cm) in diameter. The 0.25-inch (0.635 cm) section is inserted into a hole 
in the tip of the arm to attach the target to the arm. Two materials are required to create 
the target because, if only one material is used and a particle of that material impacts the 
target, the SEM would not be able to determine the particle material from the target. The 
selection of target materials to use is based on the determination of those materials that 
can withstand the plume temperature during the exposure time and that are also soft 
enough to allow sizable impact craters. Aluminum and cooper were chosen as potential 
target materials because of their softness, coupled with relatively high thermal diffusivity. 


Aluminum 



b. Front View 


Figure 5. Target Configuration 


Calculations in Chapter 3 verify that aluminum and copper can safely withstand the 
specified residence times. 



CHAPTER 3 


THERMAL RESPONSE OF THE TARGET MATERIAL 

Since the thermal conditions in a rocket exhaust plume are unfavorable for most 
common engineering materials, it is desirable to know how long the probe can reside in 
the plume before the probe experiences thermal failure. This residence time represents an 
operational limit for the PICS. 

Since the probe is composed of two materials with different thermal properties, it 
is necessary to determine the time required for each material to reach the service 
temperature (critical temperature), taken to be 75% of its melting temperature for this 
application. The smaller of these two critical times will be the design residence time in 
the plume for the probe. The calculation of the critical times will assume that the entire 
target is of a single material. It is a conservative assumption in that the calculated 
minimum critical time will be less than the value when the second material is present. 

Statement of the Thermal Problem 

The thermal response of the target as a function of time is obtained from the 
unsteady heat conduction equation which, in cylindrical coordinates, is (Myers, 1987) 

i ar a 2 T 1 aT 1 a 2 T a 2 T 

a dt dr 2 + r Sr r 2 d6 2 dz 2 
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where T = temperature 
t = time 

a = thermal conductivity 
r = distance in the radial direction 
z = distance in the axial direction 
0 = distance in the angular direction. 

Equation 1 requires two boundary conditions in each direction (z, r, and 0 ) and an 
initial condition. Figure 6 presents the target geometry. The boundary conditions for the 
target thermal problem follow: 

1. On the exposed flat face, there is a shock layer downstream of the near-normal 
portion of the bow shock. The temperature in this region is conservatively taken to be 
the static temperature downstream of the normal shock with an assumed heat transfer 
coefficient of h 2 = 1,800 B TU/hr ■ ft “ • ° R (10,220 Watt/m 2 -K), as shown in Figure 7. The 
assumed flow properties in the nozzle exhaust plume are taken from Callens (1994). 


R 


Front View Side View 




Figure 6. Target Geometry 
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The boundary condition is 


h 2 A(T 2 -T s ) = -k t A 



z = 0 , 0<r <R (2) 



2. On the curved exposed surfaces, the temperature is taken to be the value 
corresponding to the conditions after the flow expands around the comer until the 
static pressure is equal to the static pressure at the exhaust plane of the rocket. The 
ratio of heat transfer coefficients, hjfa, varies approximately with (P/P?) 1 ' 2 (Bailey, 
1965). The two boundary conditions are 

0< z< L, , r = R 


h 3 A(T 3 - T s ) = -k t A— 


R 


( 3 ) 
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h 3 A(T 3 -T,) = -k w A^ ; 


z> L, , r = R 


(4) 


In reality, the temperature is greater on the windward side of the target and probe arm 
as it travels through the plume because the pressure is greater on this side. However, 
smce the axial flow velocity is expected to be approximately 8,000 ft/sec 
(2,438 m/sec) and the target transverse velocity is expected to be approximately 
100 ft/sec (30 m/sec), this pressure difference between the windward and leeward 
sides is small. Likewise, the temperature difference is correspondingly small and 
therefore neglected. 

3. On the surfaces of the target touching the support material, the heat flow by 
conduction is assumed to be equal to the heat flow out of the surface into the arm. 
The three boundaiy conditions are 


k - A f 


, 01 . 
k ' A &' 


_ kA Si 

‘ 0Z 1 


-k A— 
sup 3z 


'Ll 


ST 

supA Sz 


«L, 


-k A— -1 
• A Sz 


r > r, , z = L, 


r < r, , z = L, 


L, <z<L 2 , r = r, 


(5) 


( 6 ) 


(7) 


4. A symmetry boundaiy condition exists along the centerline of the target and steel 
support because, in this problem, there is no temperature variation with d. The 
boundary condition along the centerline is 
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— = 0 ; r = 0 (8) 

5. The area downstream in the steel holder is considered to be insulated. Therefore, the 
boundary condition is 

dT 

— = 0 ; z = 3L 2 (9) 

6. The initial condition is that the temperature of the target and support materials is 
initially the ambient temperature which is assumed to be 540°R (300 K). 

T; = 540°R ; 0< r<R , 0<z<oo (10) 


Properties at the Nozzle Exhaust Plane 

The Mach number Mj, temperature Tj, and pressure Pj at the exhaust plane were 
calculated for an exit plane to throat area ratio of 6 (Callens, 1994) using the equations 
for isentropic flow of a perfect gas (John, 1984): 


A* 


Miff 


( 11 ) 


ul 


T = 

1 f y-l 

1 £*yr 2 


1+1 


p,= 


“Ol 


T-lwrV-' 


1 + M 


( 12 ) 


(13) 
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where A c /A* = cross-sectional area of the exit plane/cross-sectional 
area of the throat = 6 
y = ratio of specific heats 
Mi = mach number at the exhaust exit plane 
T] = temperature at the nozzle exit plane 

T 0 i = stagnation temperature at the nozzle exit plane = 6,800° R 

Pi = pressure at the nozzle exit plane 
Poi = stagnation pressure at the nozzle exit plane = 473 psia. 

The assumed values are taken from Callens (1994). The results are Mi = 2.92, 
P, = 1 1.8 psia (81 .4 kPa), and Ti = 3404°R (1 891 K). 


Properties Downstream of the Bow Shock 

The properties downstream of the bow shock are determined from the normal 
shock relations (John, 1984): 


P o2 = P o! 


r*i M , ' 
2 M ' 

H 

l + ^M, 2 
2 J 



I 7-1 
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(14) 
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T, = 


T^l + ^yV 2 

1 + ^M 2 2 


(16) 


P = P 

r 2 r l 


1 + yM , 2 
U + yM 2 : 


(17) 


The properties behind the shock wave and directly in front of the target are 
P o2 = 120.5 psia (830.8 kPa), M 2 = 0.43, P 2 = 108.1 psia (745.3 kPa), and T 2 = 6187°R 
(3437 K). 


Properties Behind the Shock Wave 
Along the Side of the Target 

It is assumed that the static pressure on the side of the target is approximately 
equal to the static pressure at the exit plane of the nozzle ( Pf=P /). The equations for 
isentropic flow of a perfect gas were used to calculate M 3 and T 3 . 




1 + ^yV 1 

u4V 


7-1 


(18) 


T 3 =T, 


f y — 1 N 


(19) 


Equation (18) was used to calculate M 3 = 2.18 and Eqn. (19) was used to calculate 
T3=4,260°R (2,366 K) with P 3 = Pj = 1 1.8 psia (81.4 kPa). Figure 8 is a summary of the 
thermodynamic properties in the flow field and around the target. 
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Semi-Infinite Anproyimafimi 

of Target Material 


Since small residence times are expected, a semi-infinite slab approximation of 
the target material was used to get an estimate of the safe residence time. This estimate 
assumes that the target and support material are made of the same material and that the 


curved surfaces are insulated . The approximation is given by Ho lman (1990): 


T - T 

s 1 

T -T 

co i 


h 2 qt 

= l-e k2 


f 


1 - erf 



( 20 ) 


where T s = critical temperature 

T; = initial temperature = 540°R 

To, = environment temperature = 6,187°R 
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h = heat transfer coefficient = 1,800 BTU/hr-ft 2 -°R 
a = thermal diffusivity of target material 
k = thermal conductivity of target material 
t = safe residence time. 

Estimates of the safe residence time were made for aluminum and copper. The 
thermal properties of these materials and the steel probe arm are listed in Table 2. 
Al umin um and copper were chosen as target materials because of their relatively high 
thermal diffusivity, low hardness, and low cost. High thermal diffusivity is required 
because of the high rate of heat transfer and high environment temperature. A high 
thermal diffusivity allows a material to dissipate heat from a point quickly. Low hardness 
is required because it is desirable to have deep, distinct craters instead of shallow craters 
that are difficult to locate. Safe residence times for aluminum and copper were estimated 
to be 77 ms and 784 ms, respectively. The importance of this calculation is that the safe 
residence time for the actual probe, with convective heat transfer to the flat and curved 
surfaces, must be less than 77 ms. 


Table 2. Properties of Materials Used in Thermal Calculations 



Density 

lbm/ft 3 

Specific Heat 
BTU/lbm°R 

Conductivity 

BTU/hrft°R 

Critical 

Temperature 

°R 

Steel 

487 

112xlQ' 3 

. 24.8 

1960 

Copper 

559 

917X1CT 4 

223 

1830 

Aluminum 

167 

207x1 O' 3 

127 

1260 
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Finite-Difference Formulation of 
Thermal Response Problem 


An exact solution of Eqn. (1) with the boundary conditions and initial condition 
given by Eqn. (2) through Eqn. (10) is difficult because of the geometry. A finite- 
difference approximation is used instead with the node geometry shown in Figure 9. 



tr 




222 




!hq 

syu 


it 


1 MMHi 




H 

3 

3 

333 

I 


Side View 


Figure 9. 


IS! /irl Ck nnAmofi- 

A vouv Nj WUtlVU 


y for t imte-Biffereacc Formulation 


The approximation uses the electrical analogy wherein thermal resistances between nodes 
and thermal capacitances of the nodal elements are used to model the flow of heat into or 
out of the node as if it were a flow of electricity as shown in Figure 1 0. 
thermal resistance to conduction and convection are modeled as 


R 


conduction 


L_ 

kA 


R 


convection 


_ 1 _ 

hA 


where L - length in the conduction direction 

k = thermal conductivity of the conducting material 
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h = heat transfer coefficient 
A = conduction or convection area. 


To Flat 

Exposed 

Surface 


To Curved Exposed Surface 



Rf. 


Tr-„ 





To Insulated 
Condition 


To Centerline of Arm and Target 

Figure 10. Electrical Analogy for Thermal Conduction Through a Node 


Thermal capacitance is modeled as 

C = pcV 

where p = density of the material 

c = specific heat of the material 
V = volume of the node. 

Nodes on a material interface require that the node be separated into a section for 
each material. Figure 11 illustrates a nodal volume where the z axis is directed 
perpendicular to the plane of the paper. The radial cross sectional areas for region A and 


region B are 
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Area A = — 
A 2 


rAr + 


At 


Area B = — 
B 2 


, Ar 
rAr-| y 


( 21 ) 


( 22 ) 



The radial cross section area for the entire volume is 

Area A +B = toAt (23) 

The finite-difference formulation for Eqn. (1), written in terms of thermal 
resistances and thermal capacitances, is (Holman, 1990) 





t r c. 


R, 


' + ^ T: P 

At 1 


j R u At 


where p = current time 

i = current node 

j = nodes around node i 


(24) 
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Rij = thermal resistance for node ij 
Q = thermal capacitance for node i 

At = time increment 

Tj P = temperature at node i for the current time step 
T; 1 ^ 1 = temperature at node z for the next time step. 

Eqn. (24) is written for each node in the system, and the equations are solved 
simultaneously for Tj P+1 . The equations for Ry and Q for each node are presented in 
Table 3. The nomenclature used in Table 3 is defined in Figure 12. 



Front View Side View 

Figure 12. Nomenclature for Finite-Difference Formulation Parameters 


The directions of r+ (radial directed towards the curved surface), r-(radial directed 
towards the centerline), z+ (axial directed towards the insulated face ), and z- (axial 
directed towards the exposed face) are illustrated in Figure 10. 
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Table 3. Thermal Resistances and Capacitances for Specified Nodes 


a. Centerline Convection Node 
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Table 3. (continued) 

c. Centerline Target-Support Node 
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Table 3. (continued) 

e. Centerline Insulated Node 



R r _ = 00 


f. Full Node on the Exposed Flat Face 
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Table 3. (continued) 

g. Full Node in the Target Material 
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Table 3. (continued) 

i. Node Along the Longitudinal Material Interface 



j. Rear Corner Node Along the Longitudinal MstcrisI Interface 
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Table 3. (continued) 

k. Full Node in the Support Material 



I 
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Table 3. (continued) 

m. Convection Corner Node 
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Table 3. (continued) 

o. Nodes on the Material Boundary on the Curved Exterior Surface 
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p. Support Node on the Curved Exterior Surface 
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Table 3. (continued) 

q. Insulated Node on the Curved Exterior Surface 
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r. Node Along a Radial Material Boundary 
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Results of the Finite-Difference Analysis 

Temperature profiles for three inches of target and support material were 
calculated. There are 10 nodes in the radial direction and 120 nodes in the axial direction. 
The distance between nodes in both the radial and axial directions is 0.025 inches 
(0.635 cm). There is convection heat transfer into both exposed surfaces as shown in 
Figure 8. A safe residence time is calculated to be 51 ms as shown in Table 4. Table 4 
also shows that the semi-infinite approximation is not acceptable for this problem 
because of the convection heat transfer into the sides. 


Table 4. Safe Residence Times 



Semi-Infinite 

Finite-Difference 

Finite-Difference 


Approximation 

Approximation 

Approximation 


(sides insulated) 

(sides insulated) 

(actual case) 


ms 

JUIO 

ms 

Copper 

784 

754 

384 

Aluminum 

77 

79 

51 


Figures 13 through 21 present the results of the calculations. The temperature 
increases faster near the comer of the front surface, as shown in Figures 13 and 14, 
because of the heating of the comer from both the front and the side. Moreover, 
al uminum heats up faster than copper because of the lower value of thermal diffusivity 
for aluminum, 1.02xl0‘ 3 ft 2 /sec (9.48xl0' 5 m 2 /sec) compared to 1.21x10° ft7sec 
(1.12X10* 4 m 2 /sec) for copper. 

Figure 15 indicates that the high temperature region has not penetrated very 
deeply into the aluminum target when the comer of the front face reaches the critical 
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temperature. Copper shows greater penetration when the critical temperature is reached, 
as seen in Figure 16, and the influence of the lower conductivity of steel is seen in the 
slope at the centerline interface between the target and support materials. 

Figure 17 indicates that the exposed curved surface of the al uminum target heats 
up most at the edge, which agrees with Figure 13. The temperature along the curved 
exposed surface decreases sharply until the target support interface is reached. This 
temperature decrease is caused by the lower heat transfer coefficient and the lower 
temperature on the side of the target. The change in slope at the material interface is 
caused by the lower conductivity of the steel. 

Figure 18 indicates that the heat penetrates deeper into the copper target than the 
aluminum target (Figure 17) along the exposed curved surface again because of the lower 
thermal diffusivity of aluminum. Figures 19 and 20 indicate that the outer edge of the 
aluminum and copper targets heat up faster than the centerline, as would be expected. 
These figures also indicate that, when the outer edge of the exposed flat face starts to fail, 
the centerline has not failed. Therefore, if the safe residence time is exceeded, useful data 
may still be gathered from the interior of the target. Figure 21 further indicates this 
possibility by showing the difference between the centerline temperature and the comer 
temperature of the exposed flat face. The finite-difference program is found in 
Appendix A, and the data represented in the Figures 13 through 21 are tabulated in 
Appendix B. 
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Figure 17. Temperature Profiles Along Exposed Curved Surface of 
Aluminum Target and Steel Holder 
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Figure 18. Temperature Profiles Along Exposed Curved Surface of 
Copper Target and Steel Holder 
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Figure 19. Axial Temperatures for Aluminum Target and 
Steel Holder at Thermal Failure (t = 51 ms) 



Figure 20. Axial Temperatures for a Copper Target and Steel Holder 
at Thermal Failure (t = 384 ms) 










CHAPTER 4 


KINETIC RESPONSE OF PARTICLES IN THE EXHAUST 


Development of Particle Velocity 
Governing Equations 

The flow near the exit plane will typically be supersonic so that a near-normal 
bow shock will be present upstream of the target as shown in Figure 22. 



The particles that impact the surface of the target will experience some deceleration in the 
shock layer between the bow shock and the surface. It is expected that the deceleration 


44 


45 


will be large for sufficiently small particles and negligible for sufficiently large particles. 

It is necessary to model the velocities of the particles in the plume both before and 
after they encounter the shock. The particle velocity immediately upstream of the shock 
serves as a boundary condition for the determination of the velocity of the particle in the 
shock layer. This boundary condition allows the calculation of the particle impact 
velocity, which is an input to the ballistic impact computer code used to calculate the 
expected impact crater depth and diameter. 

The force on the particle is given by 

F = ^C iPg A ? V, 2 (25) 


where F 

= force on the particle 

c d 

- coefficient of drag for the particle 

Pg 

= gas density 

A p 

= cross-sectional area of the particle 

v r 

= velocity of particle relative to gas. 


Newton’s second law of motion also gives the force on the particle as 

dV 

F=m ' v 'af (26) 

where m p = mass of the particle 
V p = particle velocity 
x = axial distance. 
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These two equations are combined to give the particle velocity before it enters the bow 
shock (Callens, 1994 ): 


dV p C dPg A p 


dx 


2m 


p V 


V & 

V p" 2 V s + V 


2 


where V c = gas velocity. After manipulation, Eqn. (27) becomes 


2Pp» P „ dV, 


V. 


C d P g A p p dx p 


-V 2 -2VV -V 2 =0 


g p g 


where p p = particle density 
Up = particle volume. 

By making the substitutions 

2 Pp u p 


a = 


^dPg^P 


c = 2V 


b = -1 


d = -v= 


Eqn. (28) becomes 


dV n 

aV — — ^ + bV 2 + cV n +d = 0 
p dx p p 


(27) 


(28) 


(29) 


(30) 


This equation is used to model the particle velocity in the pl um e. 


Solution Method for the 
Particle Velocity Equation 

A finite difference approximation of Eqn. (30) is obtained by replacing V P and the 
derivative by their finite difference equivalents. The result is 
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fv +V n ^| 

p»+l Pi 

fv -V ) 

P»+i Pi 

( + b 

fv +v D > 

Pi-t-l Pi 

2 

4- C 

rv D +v„) 

Pi+i Pi 

l 2 J 

Ax J 

K 2 ) 

l 2 J 


+ d = 0 (31) 


Through manipulation, Eqn. (31) becomes the quadratic equation 


f— + -lVp +( bV p +c ) V p + 

V 2 

'b 

a 

+ cV D +2d 

VAx 2 J Pi+I v Pi / p,+1 

Pi 

2 ~ 

Axj 

Pi 


= 0 


The values for a and d should also be replaced with 

2 Pp u p 


a = 


4 gi „ + p g , A 


V +v 

8i+i Si 


(32) 


(33) 


to make the finite difference approximation valid. 

Before Eqn. (32) can be solved, p g and V g must be known at node i and i+1. 
Determining these values requires the assumption of a typical nozzle shape. One-half of 
the nozzle is assumed to be of the shape shown in Figure 23. This shape is defined by the 
relations 

Ar+A* 


A c + A * I 

A c — cosl — I +• 

A A c 2 ). \x*> 


0 < — <1 
x* 


(34) 


A 

A* 


A e - — t— * 

- e <-> 


( 7t(x-xJ) A e + A * 


COS 


lsJL<;2b. 


where x = distance from the combustion chamber to the point of interest 
x e = distance from the combustion chamber to the exhaust plane 
x* = distance x from the combustion chamber to the throat 
A = cross-sectional area at x 


A c = cross-sectional area in the combustion chamber 
A e = cross-sectional area at the exhaust plane 
A* = cross-sectional area at the throat. 
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Figure 23. Assumed Nozzle Shape 


The value of A/A* is assumed to be 6 (Callens, 1994) because this value is typical for 
hybrid rocket engines. 

The nodal values for M, V gi and p s are determined from the relations (John, 1984) 


A* 

~k 





7 + 1 

2(r-0 


f\7 \ 
g 



Y-i 

2 



-i 


(35) 


( 36 ) 




49 


— = fl +"~M 2 1 ,_1 (37) 

P, 3 2 > 

where M = mach number of the gas 

y = ratio of specific heats for the gas 

a* = speed of sound at the throat 

p t = stagnation density of the gas. 

Speed of sound at the throat is calculated from (John, 1984) 

a* = ^/yRT* (38) 

where R - exhaust gas constant = 2,093 ft 2 /sec 2 -°R (350 m 2 /sec 2 -K) (Callens, 1994) 

T* = stagnation temperature at the exit plane from Eqn. (12) evaluated at M*=l . 
The particle velocity impacting the probe is given by (Callens, 1994) 

~C d SpgA p 

V p = V^e 2m > (39) 

where C d = coefficient of drag 

V p = impact velocity of particle 

Vpe = particle velocity at the nozzle exit plane (beginning of the shock layer) 

6 = shock stand-off distance. 

The shock stand-off distance is approximately (Anderson, 1989) 

6 = 0.143rexp[^) (40) 

where r = radius of the probe. Tne shock stand-off distance, o, is shown in Figure 22. 


Drag Coefficient - C 
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The computational procedure which solves Eqn. (32) through Eqn. (40) to 
determine the mach number, gas velocity, gas density, particle velocity at the exit plane, 
and particle impact velocity at the target surface is as follows: 

1. Divide the flow region from the combustion chamber to the exhaust plane into N 
regions Ax apart. 

2. For each x i+1 , calculate A/A* from Eqn. (34), starting with known values at the 
combustion chamber. 

3. Calculate M i+ j corresponding to the area ratios from Eqn. (35). 

4. Calculate (V J i+I corresponding to the mach numbers from Eqn. (36). 

5. Calculate (Pg)t+i from Eqn. (37). 

6. For an assumed particle size, calculate (V p ) i+1 from Eqn. (32), using the above values 
of V g axid p g . Ch is given in Figure 24 (White, 1994). 

7. Once the exit velocity is known, calculate impact velocity from Eqn. (39). 



Figure 24. Drag Coefficient of a Smooth Sphere 
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Results of the Analysis 

Figures 25 through 27 show the mach number, gas velocity, and gas density, 
respectively, for the assumed nozzle shape. These graphs are valid for any length nozzle 
of the assumed shape as long as A/A* = 6 and x*/x e = 1/3. Figures 28 through 32 show 
the velocities of various sized particles of hydroxl-terminated polybutadiene 
(p p = 115.5 lbm/ft 3 , 1850 kg/m 3 ) exiting nozzles of lengths of 1, 6, 12, 24, and 48 inches 
(2.54, 15.24, 30.48, 60.96, and 121.9 cm). This material was chosen for the particles 
because it is the least dense material likely to be present in solid form in the exhaust. 


Therefore, if a particle of this material will create a usable crater, so will a particle of a 
denser material. The velocities at the exhaust plane are the only values of interest. It is 
evident from Figure 28 that particles larger than approximately 500 pm may not have 


enough momentum exiting a 1-inch (2.54 cm) nozzle to create a useiul crater, evaluation 
of the craters formed by these particles will be required to determine if this situation 
constitutes an operational limit of the PICS. All other particles on the various graphs in 
Figures 29 through 32 appear to have sufficient velocity entering the shock to form usable 
craters. 


Figure 33 shows the approximate impact velocity of particles up to 3,000 pm in 
diameter exiting from nozzles 1 to 48 inches in length. It is evident that particles less 
than approximately 0.1 pm in diameter will not make usable craters because of 
deceleration in the shock layer. Impact velocities are low for particles larger than 
3,000 pm because of lower acceleration in the nozzle. It is concluded that the PICS has 


the capability of characterizing particles between 0.1 pm and 3,000 pm in diameter ii 
nozzles greater than 1 inch (2.54 cm) in length. The particle dynamics program is ir 
Appendix C, and the data represented in the graphs are found in Appendix D. 
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Figure 25. Mach Number in the Nozzle 
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Figure 30. Particle Velocity in the Nozzle for X e = 12 inch 
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Figure 31. Particle Velocity in the Nozzle for X e = 24 inch 
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Estimation of Crater Depth and Diameter 

Crater depths and diameters, as estimated by McMurtry’s projectile penetration 
computer code (McMurtry, 1997), for various HTPB particle sizes are given in Figures 34 
through 37. This application of McMurty’s code assumes that the particles are spherical. 
Small particles are roughly spherical in the plume because the flow field around these 
particles, coupled with the high heat in the plume, tends to mold the particles into a 
sphere. As particle size increases, this spherical assumption is less accurate. Fortunately, 
only relatively small particles are expected to be in the exhaust plume. Table 5 gives the 
penetrator and target material data used for the crater analysis (Beer and Johnston, 1992). 
The HTPB data is an estimate of the properties of a hard rubber. Impact velocities for the 
various size particles are given in Figure 33. 


Table 5. Properties of Materials for Penetration Code 



HTPB 

(penetrator) 

6061-T651 

Aluminum 

(target) 

Copper 99.9% 
(target) 

Density (lbm/ft 3 ) 

74.9 

169 

537 

Compressive Strength (ksi) 

7.3 

45 

56.6 

Lateral Speed of Sound (ft/sec) 

8,169 

17,480 

12,927 

Transverse Speed of Sound (ft/sec) 

4,921 

10,007 

9,268 

Poisson’s Ratio 

0.35 

0.33 

0.33 

Hugoniot Constant 

1.577 

1.338 

1.489 
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Figures 34 through 37 show that the diameter of the crater formed in 6061-T651 
aluminum by a particular size HTPB particle is unique for a given nozzle length. 
However, error introduced by the measurement process requires the measurement of both 
the crater depth and diameter. The particle velocity then can be determined from graphs 
similar to Figures 28 through 32 for the known nozzle configuration. Longer nozzles 
allow particles to accelerate more before exiting the engine, enabling particles to create 
larger craters as shown in Figures 34 through 37. 

Figures 35 and 37 indicate that some particles in the shorter nozzles do not 
generate sufficient velocity to form a crater. This condition creates an operational limit, 
which is a function of nozzle length, on the size of the particles which can be 
characterized by the PICS. Figures 36 and 37 indicate the same trends for impact with 
copper that are seen in irnpact with slun iinuni (h'i^ures 3 A- and. 35}. Crater depth and 
diameter are slightly higher in aluminum than they are in copper because of the higher 
compressive strength of copper. 
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Figure 36. Crater Diameter in Copper as a Function of HTPB Particle Size 
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CHAPTER 5 


DESIGN OF PROTOTYPE 
Energy Analysis of the Probe 

For the design residence time in the plume of 20 ms, a target velocity of 
approximately 100 ft/sec (30.5 m/sec) is required. The target launch velocity should be 
higher than this to correct for aerodynamic drag and pressure resistance upon entrance to 
the plume. To determine the correct arm rotational launch velocity, it is required to 
calculate the kin etic energy of the arm Ej, the energy of the arm lost in the entrance to the 
pl um e E 2 , and the energy of the arm lost to drag inside the plume E 3 for a target 
translational velocity V T = r T co where r T is the radius of curvature for the path of the target 

and co is the arm rotational velocity. 

The kinetic energy of the arm is 

E, = ^1®’ (41) 

where I = mass moment of inertia of the arm 
co — angular velocity of the arm. 

The energy change in die entrance to the plume is given by 

E ,=:>,- p > (t > ds <42) 


63 


64 


where R radius of the probe at the entrance to the plume 
Pgas = gas pressure in the plume 
Patm = atmospheric pressure 
A(t) = probe profile area as a function of time 
V t = probe translational velocity 
dS = V t dt. 

In Eqn. (42) P gas - P atm and V t are both constant. Figure 38 illustrates the probe entering 
the plume. An expression for x(t) is required to determine A(t). 



The distance from the probe centerline to the plume, y(t), is y(t) - R-V t t. From the 
Pythagorean theorem, 


which reduces to 


x(t) = y'R 2 -(R- V t t) 2 


(43) 
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x(t) = A /2V t Rt-V t 2 t 2 (44) 

Substituting A (t) = 2 x(t) L p , where L p is the length of probe entering the plume, 
into Eqn. (42) yields 

2R 

E 2 - 2L„V, (? m - P„ ) ^V.Rt-V.Vdt (45) 

which reduces to 

E2 71 Lp R (Pgas " Patm) (46) 

The probe actually gains energy entering the plume because the pressure in the plume is 
1 1.8 psia (8.1 kPa) which is less than the atmospheric pressure of 14.7 psia (10.1 kPa). 

Assuming a plume width of 6 inches (15.24 cm), the energy lost to drag inside the 
plume is 

E, = (47) 

where Cd = coefficient of drag of the arm (conservatively assumed to be 2) 

Pgas = gas density in the plume 
Aprofiie = probe profile area = 2RLp 
V t = probe translational velocity. 

The launch energy is given by 

|lo L 2 =|la) p I +E 2 +E 3 ' (48) 

where cq l = rotational velocity of the arm at launch 

c o p = rotational velocity of the arm in the plume. 


66 


The last two terms are highly sensitive to the orientation of the probe arm when it 
enters the plume. If the arm enters the plume at a negative angle of attack, it will 
experience much higher energy losses than this analysis assumes. The energy lost to drag 
E 2 is also highly sensitive to whether the nozzle is underexpanded or overexpanded. 
These nozzles produce regions of high and low pressure in the plume which will affect 
the drag on the probe. Calibration of the PICS is required to correct for energy loss 
caused by these two problems. 

The mass moment of inertia of the arm is required for the calculation of the 
launch rotational velocity of the arm. The mass moment of inertia of an object about an 
axis is defined as 



(49) 


To evaluate this equation, the arm is divided into several sections as shown in Figure 39. 


Figure 39, Sections of Arm Used in Mass Moment of Inertia Calculation 

The mass moment of inertia of each section is found, and these moments are summed. 

The moments of inertia are f 

I, = {'’ S rpu(5) J dr 

I, = £r ! p<5) 2 cos(45»)-'di 

I 
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1 3 = J 3 r 2 p7r(.325) 2 cos(45°) ! dr (50) 

1 4 = ^r 2 pn(25) 2 cos(45°) _1 dr 

1 5 = -~pTc(.25) 2 (35)(25) 2 + P 7x(^25) 2 (3_5)(2 1) 2 

where p - density of steel 

r = distance from the axis of rotation. 

Summing these moments of inertia yields 527 lbm-in 2 (1.52 kg-m 3 /sec 2 ). 

The launch energy is calculated to be 185.8 lbf-ft (251.9 kg-m 2 /sec 2 ). This energy 
corresponds to an arm rotational velocity of 57.1 rad/sec and a target translational launch 
velocity of 101.3 ft/sec (30.9 m/sec). 

A similar analysis can be made for the minimum safe residence time of 50 ms and 
is calculated to be 29.6 lbf-ft (40.1 kg-m 2 /sec 2 ). This energy corresponds to an arm 
rotational velocity of 22.8 rad/sec and a target translational launch velocity of 39.9 ft/sec 
(12.2 m/sec). The results of this analysis for the design condition are listed in Table 6. It 
is observed that the calculated losses are negligible for the assumed conditions. 

Design of the Spring 

The spring is viewed as a deflected beam where the energy in is given by (Juvinall 
and Marshek, 1991): 

f i T 2 

Energy = J s 


( 51 ) 
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where L = length of the beam 
E = Young’s modulus 
I = area moment of inertia 
T = torque applied to the beam 
x = distance along the centerline of the beam. 


Table 6. Velocities and Energies the Probe Encounters During Operation 


Rotational Velocity at 

Launch 57 rad/sec 

©L 

Energy of the arm at 

Launch 185.8 lbf-ft 

E l 

Rotational Velocity in 

the Plume 57 rad/sec 

© D 

Energy Change Entering 

the Plume 0.17 lbf-ft 

Epe 

Translational Velocity 

of Target Launch 99.9 ft/sec 

v L 

Energy Change to Drag 

for the Entire Path -0.033 lbf-ft 

Epd 

Tm hhIa^iam A i X 7 /\i jr\ nifir 
i i auMauuum r ciuwuj 

of T arget in the 1 00 ft/sec 

Plume 

Ve 

TT n nrmr At tn o 
XVllI VUW JUUVt VA WV 

Arm in the Plume 185.7 lbf-ft 

E P 


Ebsco Spring Company, Inc., produces software which calculates, among other 
things, the torque a torsion spring will produce as a function of the angular displacement 
of the spring arm. It also calculates the safe stress, which includes an empirically 
determined safety factor, for a particular spring and the stress in the spring as a function 
of angular displacement. The program is called Spring Calculation Program and is 
available on the Internet at “http://www.tulsaweb.com/ebsco.” Inputs and outputs, used 
in the spring design, for Spring Calculation Program are given in Table 7. 
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Table 7. Inputs and Outputs for Spring Calculation Program 


Inputs 

Outputs 

1 . Spring Material (The program lists 

1 . Unbended Length, L in Eqn. (5 1 ), 

material choices and their relevant 

of the Spring Wire 

properties) 


2. Mean Diameter of Coils 

2. Torque, T in Eqn. (5 1), in the Spring 


for Each of the Deflection Angles 

3. Leg Length 

3. Stress in the Spring for Each of the 


Deflection Angles 

4. Number of Coils 

4. Safe Working Stress 


in the Spring 

5. Spring Deflection Angles 



Iteration between Eqn. (51) and the Spring Design Program is required to design 
the torsion spring. A spring is acceptable if the unbended length L of spring wire and the 
torque T produced by the spring satisfies Eqn. (51) when Energy = 186 Ibfft 
(252 kg-m 2 /sec 2 ). The design torque is calculated to be 2275 lbf-in (169.5 kg-m 2 /sec 2 ) 
and the unbended length of wire is 195 inches (495.3 cm). This spring fails when the 
stress in the spring exceeds 176,250 psi (1.215 GPa). The specified spring is given in 
Table 8. 

Figure 40 indicates the torque of the specified spring versus angular deflection of 
the spring arm. The design angular deflection of the spring arm to produce the required 
20 ms residence time is approximately 106°. Figure 41 indicates that the spring can 
theoretically be calibrated to produce target launch velocities between 0 and 140 ft/sec 
(61.0 m/sec) by varying the angular deflection of the spring arm. It is evident from 


/ 
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Figure 42 that no test requirement will require a deflection which will exceed the safe 
working stress of the spring. 


Table 8. Specifications for Required Torsion Spring 


Number of Coils 

11.125 

Wire Size 

0.625 inch 

Mean Diameter 

5 inch 

Material 

Chrome Silicon 

Leg Length 

6 inch 

Body Length 

Min 7.813 inch 


Max 8.273 inch 

Max Shaft Size 

4.021 inch 

Safe Working Stress 

176,250 psi 



Figure 40. Spring Torque Versus Spring Arm Deflection 
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Figure 41. Target Translational and Rotational Velocity 
Versus Spring Arm Deflection 



Figure 42. Stress in Spring Versus Spring Arm Deflection 
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Parts List and Shop Drawings 

Table 9 is a list and description of parts required to build the PICS. 


Table 9. Parts List 


Figure 


Description of Part 


Figures 43 through 52 are Required to Manufacture 
the Probe and Spring Structure 


43 

44 

45 


Probe Arm 

Figure 43 is a drawing of the unbended length of rod required to produce the 
probe arm, and Figure 44 shows the required bends to produce the final shape 
of the probe arm. The 0.25-inch (0.635 cm) hole which holds the target is 
shown in Figure 45. ASTM A-514 steel is specified for this part because of 
its hardness. 


Probe Collar 

46 The probe collar attaches the probe arm to the shaft. It is suggested that the 
vertical and horizontal sections of the probe collar be manufactured separately 
and welded together; Weld size is not critical on this component since it is a 
pinned joint and has little mass moment of inertia of its own. ASTM A-36 
steel is specified for this part. The probe arm should be inserted into the 
probe collar and welded in place. This was not done on the prototype because 
the probe collar was designed after the probe arm was built. It was later 
realized that welding the parts would create a mechanism with less play than a 
threaded arm and probe collar would create. 

Shaft 

47 The shaft holds the spring mandrel and probe collar in place. Nuts axe 
required on this part only to keep the shaft from slipping out of the vertical 
supports. Material is ASTM A-36 steel. 
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Vertical Support 

Vertical supports hold the shaft in place. The distance between the bottom of 
the vertical support and the hole is critical. This distance is chosen because it 
allows the spring to sit on the mandrel properly. Material is AS nvl A-36 steel. 
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Baseplate 

The baseplate holds the vertical supports, horizontal supports, and spring arm 
holder in place. Material is ASTM A-36 steel. 
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Table 9. 
Figure 

50 

51 


(continued) 

Description of Part 

Mandrel 

The mandrel provides support for the spring while it is deflected. A 4-inch 
(10.16 cm) outer diameter steel pipe may be used to manufacture the mandrel 
shown in Figure 50. The mandrel end plugs, Figure 51, are inserted into the 
ends of the mandrel and spot welded in place. They allow the mandrel to sit 
properly on the shaft. Material is ASTM A-36 steel. 


52 


53 


54 


55 


56 


57 
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Spring Arm Holder 

The spring arm holder keeps the spring from rotating around the mandrel. 
Material is ASTM A-36 steel. 

Figures 53 through 61 are Required to Manufacture the Pneumatic Clip 

Horizontal Support 

Horizontal supports hold the pneumatic clip in place. Material is ASTM A-36 
steel. 

Slider Holder 

The slide holder creates a pocket which the slider slips into before the clip is 
pressurized. Material is ASTM A-36 steel. 

Pin Holder 

The pin holder holds the pin in place above the slider. Material is ASTM 
A-36 steel. 


Pin 

The pin keeps the slider from slipping out of the slider holder while the 
pneumatic clip is pressurized. Material is ASTM A-5 14 steel. 

Bottom Plate 

The bottom plate holds the pin, pin holder, and slider holder in place. A 
0.0625-inch (15.8 mm) rubber gasket is placed in the outer groove. Material is 
ASTM A-36 steel. 

Top Plate 

The top plate holds the pressure inlet and electric valve in place. A 
0.0625-inch (15.8 mm) rubber gasket is placed in the outer groove. Material is 
ASTM A-36 steel. 
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Table 9. (continued) 
Figure 

59 The pressure 
A-36 steel. 


Description of Part 

Pressure Vessel 

vessel is made of 4-inch (10.16 cm) pipe. Material is A STM 
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Plunger 

The plunger pushes the pin down so that it interferes with the slider and holds 
it in place while the pneumatic clip is pressurized. The two pieces should be 
made separately and welded together. A 0.0625-inch (15.8 mm) rubber 
gasket is placed on the outer race between the plunger and bottom plate. 
Material is ASTM A-36 steel. 
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Slider 

The slider holds the arm back when the pneumatic clip is pressurized. 
Material is ASTM A-36 steel. 


Miscellaneous Parts Which Do Not Require Part Drawings 


torsion spring with specifications listed in die previous section (one part) 

0.3 75-inch (0.953 cm) diameter bolts, each 4 inches (10.16 cm) in length (five parts) 
0.375-inch nuts (five parts) 

A vtrool'itfvro Afrvr** -niQTTcA 
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0.25-inch (0.635 cm) diameter bolt, each 2 inches (5.08 cm) in length (one part) 
0.25-inch nut (one part) 


3- inch (7.62 cm) spring clip (one part) 

4- inch (10.16 cm) quick link (one part) 

24 inches (60.96 cm) of 0.375-inch (0.953 cm) chain 

1.5-inch (3.81 cm) width x 3.5-inch (8.89 cm) height x 0.12-inch (0.3175 cm) wire 
diameter compression spring (one part) 


0.125 inch (0.3175 cm) threaded air inlet (one part) 

0.125-inch nipple (one part) 

12 V DC, 0.125-inch connector, electric valve (75 psig = 517 kPa max) (one part) 

0.125-inch (0.3 175 cm) thick rubber gasket, 

ID: 4.00 inches (10.16 cm) OD: 4.25 inches (10.795 cm) (one part) 

0. 125-inch thick rubber gasket, 

ID: 3.50 inches (8.89 cm) OD: 4.00 inches (10.16 cm) (one part) 
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Figures 62 through 64 indicate which parts must be welded together to complete 
the fabrications of the individual pieces of the PICS. Figure 62 shows how to weld the 
slider holder, pin holder, and bottom plate together. The pin must be placed inside the 
pin holder before the parts are welded together. Figures 63 and 64 indicate how to weld 
the horizontal and vertical supports, and the spring arm holder to the baseplate. 


77 





I 

K 

1 



Figure 43. Section Required to Manufacture Probe Arm 


3.375 
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Figure 45. Call-Out of Hole in Tip of Probe Arm 
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Figure 47. Shaft 




7.50 + 0.00 



Figure 49. Baseplate 












Figure 54. Slider Holder 




AH dimensions in inches 
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Figure 55. Pin Holder 




Unless otherwise noted all tolerances 



Figure 56; Pin 





Figure 58 . Top Plate 




Figure 59. Pressure Vessel 
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Figure 61. Slider 
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Figure 64. Welds Required to Attach Vertieal Supports and Spring Arm Holder Baseplate (Front View) 
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Prototype Assembly 

Once all parts are welded together, they must be assembled. The first step in 
assembly is to place the mandrel into the spring as shown in Figure 65. Then, insert the 
shaft through the left vertical support, through the probe collar, and into the middle 
vertical support (Figure 66). Next, insert the straight spring arm into the spring arm 
holder (Figure 67), and push the spring into place. Then, insert the shaft through the 
mandrel and into the right vertical support. 



Figure 65. Spring with Mandrel Partially Inserted 


The pneumatic clip must now be assembled. First, place the 4-inch OD gaskets 
into the outer grooves in the top and bottom plates (Figure 68). Then, place the pressure 
vessel into the outer groove of the bottom plate. Place the re mainin g gasket into the 
pressure vessel, followed by the plunger and the compression spring (Figure 69). The top 
plate now can be attached to the bottom plate by means of the 0.375-inch (0.9525 cm) 
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bolts, washers, and nuts (Figure 70). Tighten the nuts, making sure that the pressure 
vessel mates with the outer groove of the top plate. 

The pneumatic clip is now ready to be attached to the baseplate. Insert the 
remaining bolt through one end of the chain Next, insert the bolt through a link on the 
chain, through one hole on the horizontal support, through the hole in the pin holder, and 
into the hole in the other horizontal support (Figure 71). Put a nut on the other end of this 
bolt. Put the quicklink onto the other end of the chain. The PICS is now assembled. 



Figure 66. Shaft Inserted Through the Probe Arm Collar 
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Figure 71. Pneumatic Clip Attached to Baseplate 


Prototype Testing 

The prototype was tested in the Louisiana Tech University Mechanical Laboratory 


to determine the velocity or the probe arm as a function of angular displacement of the 


spring arm from the lannch nrtotfirm fn nlnma ' 1 U « i • ,i i • « . . 
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chronograph designed to measure the velocity of a bullet which passes over two sensors. 
When the bullet passes over the first sensor, the chronograph activates; as the bullet 
passes over the second sensor, the chronograph deactivates. The chronograph divides the 
distance between the sensors by the time of activation to determine the velocity of the 
bullet. Since the probe arm has an angular velocity which consists of both x and y 
components, the velocity was measured while the probe arm was vertical (Figure 72). 
This configuration eliminated the y-component of velocity and resulted in a direct reading 
of the translational velocity of the probe arm. This configuration also limited measurable 


data to between the angles of 0° and 90°. The sensors were separated by 1.5 inches 
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This moment provides extra energy to the system. However, at full deflection this 
moment adds only 7 lbf-ft of energy to the system. 



Figure 73. Experimental and Theoretical Probe Velocity 


Operation of the PICS 

1. Before the test firing of the rocket, mount the PICS in a position which allows the 
probe arm to pass perpendicularly to the exhaust flow'. The probe should not be in the 
plume before launch nor should it be in the plume after launch, and the hole in the 
0.5-inch (1.27 cm) section of the probe arm should be facing the rocket exhaust. The 
PICS should be firmly secured to the substrate with either bolts or clamps. 

2. Push the slider into the pneumatic clip until it clicks into place. 

3. Pressurize the pneumatic clip to 70 psig (482.6 kPa). 
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with the jig. At this point, the 0.5-inch diameter section will be a half-cylinder. 

6. Remove jig part 2, place jig part 3 (Figure 77) over the flat surfaces of the target, and 
secure it to jig part 1 . 

7. At this point, the 0.25-inch (0.635 cm) diameter should be exposed in the center on 
the jig. Use an endmill to remove the exposed areas of the target so that the target is 
flush with the jig. 

8. The aluminum half of the target is now fabricated. Remove it from the jig. 

9. Repeat steps 1-8 with a copper (pure) cylinder. 

The jig manufactured in the Louisiana Tech machine shop and shown in Figures 

75 through 77 allows for up to four target halves to be manufactured at one time. It is a 


simple matter to increase this number by making the jig longer and adding more cradles 

fnr taraptc if* Rapquoa nf* A 4 — 
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and copper, it is suggested that only one material type be placed in the jig at one time. 


I 


fi 
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figure 74. Insert Required to Manufacture Targets 




0.875 
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/* This program calculates the safe residence time of the target and probe. If 
SEMI_INFINITE is defined the calculation assumes that the sides axe insulated. Either 
AL (aluminum) or CU (copper) must be defined for the program to work correctly.*/ 

#include <stdio.h> 

#include <conio.h> 

^include <math.h> 

#include <stdlib.h> 

#include <string.h> 

#defme INFINITE 1 .0e+300 

#define NUMR 1 1 
#defineNUMZ 121 
#define EPS 0.01 

#define dSEMI_INFINITE 
#define CU 

double T[NUMR][NUMZ]Jold[NUMR][NUMZ],To; 
double Rip,Rim,Rjp,Rjm,SumR,8umTR; 
double dt,dr,dz,t,C; 
int ij, converged; 

void main(void) { 

double hl,h2,Tl,T2; 
int pl,p2,p3,l; 
double cs,denss,ks,Tcrits; 
double ct,denst,kt,Tcritt; 
double dq.dv; 
double pi,q,r; 

char target[12],FrontTemp[20],SideTemp[20],CenterTemp[20],Dim[3]; 

FILE *Side, *Center, *Front; 

void CheckF orCritical(doub!e, double, int, int); 

void GaleulateNode(double, double, double, double); 

clrscr(); 

_setcursortype(_NOCURS OR) ; 

pi=4*atan(l); 

t=0; 


dt=.001; 

/* s 

*/ 

dr=.00635/(NUMR-l); 

/* m 

*/ 

dz=.0762/(NUMZ-l); 

/* m 

*/ 


pl=(NUMZ-l)/6.0; 

p2=(NUMR-l)/2.0; 

p3=(NUMZ-l)/3.0; 


/*** load environment data ***/ 

Tl=3437; /* K *1 

T2=2224; /* K */ 

hi =10220; /* W/m A 2*°C */ 

h2=3371; /* W/m A 2*°C */ 

/*** load material data ***/ 

#ifdef AL 
strcpy(target, "Al"); 
strcpy(FrontTemp, "alfront.dat"); 
strcpy(SideTemp, "alside.dat"); 
strcpy(CenterTemp, "alcenter.dat"); 
denst=2700; /* kg/m A 3 */ 

ct=866; /* J/kg*°C */ 

kt=220; /* W/m*°C */ 

Tcritt=700; /* K */ 

#endif 


#ifdefCU 
strepy(target,"Cu"); 
strcpy(FrontTemp,"cufront.dat"); 
strc Py(SideTemp, "cuside.dat"); 
strcpy(CenterTemp, "cucenter.dat"); 


denst=8954; 

ct=384; 

kt=386; 

Tcritt=1017; 

#endif 


*/ 

*/ 


/* kg/m A 3 
/* J/kg*°C 
/* W/m*°C */ 
/* K V 


#ifdef STEEL 

strcpy(target, "Steel "); 

denst=7800; 

ct=470; 

kt=43; 

Tcritt=1360; 

#endif 


/* kg/m A 3 */ 
/* J/kg*°C */ 
/* W/m*°C */ 
/* K */. 


denss=7800; /* kg/m A 3 */ 

cs=470; /* J/kg*°C */ 
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ks=43; /* W/m*°C */ 

Tcrits=1360; /* K */ 

Front=fopen(FrontTemp,"wt"); 

Side=fopen(SideTemp,"wt"); 

Center=fopen(CenterT emp, "wt"); 

fprintf(Side," "); 
f^rintf(Center," "); 
f^rintf(Front," "); 
for(i=0;i<NUMZ;i++) { 
fprintf(Side,"%4.21f ",i*dz*39.37); 

§>rintf(Center, "%4.21f ",i*dz*3 9.37); 

} 

for(j=0;j<NUMRy++) fprintf(Front,"%4.21f "j*dr*39.37); 

strcpy(Dim, "2D "); 

#ifdef SEM JNEINITE 
strcpy(Dim,"lD"); 

#endif 

printf "target material: %s - %s conduction\n", target, Dim); 
printf("dr=%1.2e in, dz=%1.2e in, dt=%1.2e sec\n",dr*39.37,dz*39.37,dt); 

for(i=0;i<NUMR;i++) 

for(j=Oy<NUMZg++){ 

T[i]0]=300; /* K */ 

Told[i][j]=T[i]0]; 

} 

do{ 

t+=dt; 

do{ 

converged=l; 
for(j=NXJMZ-l y>=Oy— ) 
for(i=NUMR- 1 ;i>=0;i--) { 
r=i*dr; 

To=T[i3jj]; 

if(i— 0){ /* along the centerline */ 

if(j=0){ /* flat exposed surface */ 

Rip=dr/((r+dr/2)*pi*(dz/2)*kt); 

Rim=INFINITE; 

Rjm=l/((pi/8)*dr*dr*hl); 
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Rjp=dz/((pi/8)*dr*dr*kt); 

dv=(pi/8)*dr*dr*(dz/2); 

C=denst*et*dv; 

_ CalculateNode(T[i+l][j],T[i+l]jj],T[i]0+l],Tl); 
}else if(j— NUMZ-1){ I* insulated surface */ 

Rip=dr/((r+dr/2)*pi*(dz/2)*ks); 

Rim=INFINITE; 

Rjm=(dz)/((pi/8)*dr*dr*ks); 

Rjp=INFINITE; 

dv=(pi/8)*dr*dr*(dz/2); 

C=denss*cs*dv; 

CalculateNode(T[i+l]0] 5 T[i+l]0] 5 0,T[i][j-l]); 

}else ifij=p3){ /* comer between target and arm */ 

Rip=dr/((r+dr/2) * dz/2) * ( 1 /ks+ 1 /kt) ; 

Rim=lNFINITE; 

Rjm=dz/((pi/8)*dr*dr*kt); 

Rjp=dz/((pi/8)*dr*dr*ks); 

dv=(pi/8)*dr*dr*dz/2; 

C=(denss*cs + denst*ct)*dv; 

CalculateNode(T[i+ 1 ] jj] ,T[i+l ] [j],T[i] [j+ 1 ],T[i] 0- 1 ]); 
}else if(j<p3){ /* interior target nodes *1 

Rip=dr/((i+dr/2)*pi*dz*kt); 

Rim=INFINITE; 

Rjm=dz/((pi/8)*dr*dr*kt); 

Rjp=dz/((pi/8)*dr*dr*kt); 

dv=(pi/8)*dr*dr*dz; 

C=denst*ct*dv; 

CalculateNode(T[i+ 1 ] 0] ,T[i+l ] 03 >T[i] [j+ 1 ] ,T[i] [j- 1 ]); 
}eise if(j>p3){ /* interior steel nodes */ 

Rip=dr/ ((r+dr/2)*pi*dz/2 *ks); 

Rim=3NFINITE; 

Rjm=dz / ((pi/8)*dr*dr*ks); 

Rjp=dz/((pi/8)*dr*dr*ks); 

dv=(pi/8)*dr*dr*dz; 

C=denss*cs*dv; 

^ CalculateNode(T[i+l ] [j],T[i+ 1 ] U],T[i] [j+ 1 ] ,T[i] jj- 1 ]); 

}else if(j=0){ /* along flat face */ 

if(i=NUMR- 1 ) { /* comer node */ 

Rip=l/(pi*r*(dz/2)*h2); 

Rim=dr/ ((r-dr/2)*pi * (dz/2) *kt) ; 

#ifdef SEMI_INFINITE 
Rip=u\TuNTTE; 


#endif 

Rjm=l/(.5*pi*(r* dr-dr* dr/4)*hl); 

Rjp=(dz)/i (. 5 *pi *(r * dr-dr * dr/4) * kt) ; 
dv=.5*pi*(r*dr-dr*dr/4)*(dz/2); 

C=denst*ct*dv; 

CalculateNode(T2,T[i- 1 ] [j],T[i][j+l],Tl); 

}else{ /* exterior convection nodes */ 

Rip=dr/((r+-dr/2)*pi*(dz/2)*kt); 
Rim=dr/((r-dr/2)*pi*(dz/2)*kt); 
Rjm=l/(r*dr*pi*hl); 

Rjp=(dz)/(r*dr*pi*kt); 

dv=r*dr*pi*(dz/2); 

C=denst*ct*dv; 

CalculateNode(T[i+l]D],T[i-l]D],T[i]D+l] 5 Tl); 

} 

}eise if(i=NUMR-l){ /* along curved exposed surface */ 
if(j==NUMZ-l){ /* insulated comer */ 

Rip= 1 /(pi*r*(dz/2)*h2); 
Rim=dr/((r-dr/2)*pi*(dz/2)*ks); 

#ifdef SEMI_INFINITE 
Rip=INFINITE; 

#endif 

Rjm=dz/(.5 *pi* (r* dr-dr * dr/4) * ks) ; 
Rjp=INFINITE; 

dv=.5 *pi*(r*dr-dr*dr/4)*(dz/2); 

C=denss*cs*dv; 

CalculateNode(T2,T[i- 1 ] D],0,T[i] [j- 1 ]); 

}else if(j=pl){ /* boundary between materials */ 

Rip= 1 /(pi*r*dz*h2); 

Rim=dr/(((r-dr/2)*pi*(dz/2))*(l /ks+ 1/kt)); 

#ifdef SEMI_INPINITE 
Rip=INFINITE; 

#endif 

Rjm=dz/(.5*pi*(r*dr-dr*dr/4)*kt); 
Rjp=dz/(.5*pi*(r*dr-dr*dr/4)*ks); 
dv=.5*pi*(r* dr-dr* dr/4)*dz; 

C=.5*dens$*cs*dv + ,5*denst*ct*dv; 
CalculateNode(T2,T[i- 1 ] {j],T[i] [j+ 1 ],T[i] D- 1 ]); 

} else if(j<p 1 ) { /* target convection nodes */ 

Rip=l/(pi*r*dz*h2); 

Rim=dr/((r-dr/2)*pi*dz*kt); 

#ifdef SEMI_INFINITE 
Rip=INFINITE; 
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#endif 

Rjm=dz/(.5 *pi*(r*dr-dr*dr/4)*kt); 
Rjp=dz/(.5*pi*(r*dr-dr*dr/4)*kt); 
dv= 5*pi*(r*dr-dr*dr/4)*dz; 

C=denst*ct*dv; 

CalculateNode(T2,T[i-l][j],T[i][j+l],T[i][j-l]); 

}else if(j>pl){ /* steel convection nodes */ 

Rip=l/(pi*r*dz*h2); 

Rim=dr/((r-dr/2)*pi*dz*ks); 

#ifdef SEMI_INFINITE 
Rip=INFINITE; 

#endif 

Rjm=dz/(.5*pi*(r*dr-dr*dr/4)*ks); 

Rjp=dz/ (. 5 *pi *(r* dr-dr* dr/4) *ks); 
dv=.5*pi*(r*dr-dr*dr/4)*dz; 

C=denss*cs*dv; 

^ CalculateNode(T2,T[i-l]0],T[i][j+l],T[i][j-l]); 

} else if(i==p2) { /* along longitudinal boundary */ 
if(j=pl){ /* front comer of stem */ 

Rip=dr/ ((r+dr/2)*pi*dz/2)*( 1 /ks+ 1 /kt); 

Rim=dr/ ((r-dr/2) *pi *dz*kt); 

Rjm=dz/(r*dr*pi*kt); 

Rjp=dz/ ((pi/2 *(dr * dr/4+r *dr) *ks)+(pi/2 * (r * dr- 
dr*dr/4)*kt)); 

C=denst*ct*pi*r*dr*(dz/2) + 

denss*cs*pi/2*(r*dr+dr*dr/4)*dz/2 + 
denst*ct*pi/2*(r* dr-dr* dr/4)*dz/2; 

CalculateNode(T[i-fl]0],T[i-l]D],T[i]0+l],T[i]D-l]); 

}else if(j=p3){ /* rear comer of stem */ 

Rip=dr/ ((r+dr/2) *pi*dz*ks); 

Rim=dr/ ((r-dr/2) *pi *dz/2) *( 1 /ks+ 1 /kt); 

Rjm=dz/((pi/2*(dr*dr/4+r*dr)*ks)+(pi/2*(r*dr- 

dr*dr/4)*kt)); 

Rjp=dz/(r*dr*pi*ks); 
C=denss*cs*pi/2*(r*dr+dr*dr/4)*dz + 

.5*denss*cs*pi/2*(r* dr-dr* dr/4)*dz + 
.5*denst*ct*ni/2*(r* dr-dr* dr/4)*dz; 
CalculateNode(T[i+l]0] 5 T[i-l][j] s T[i]0+l],T[i][j-l]); 

}else if(j<pl){ /* target nodes */ 

Rip=dr/ ((r+dr/2)*pi*dz*kt); 
Rim=dr/((r-dr/2)*pi*dz*kt); 

Rjm=dz/i (r*dr*pi*kt); 
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Rjp=dz/(r*dr*pi*kt); 

dv=T*dr*pi*dz; 

C=denst*ct*dv; 

CalcuiateNode(T[i+l][j],T[i-l][j],T[i]0+l],T[i]0-l]); 

}else if(j==NUMZ-l){ /* insulated face */ 

Rip=dr/((r+dr/2)*pi*dz/2*ks); 
Rim=dr/((r-dr/2)*pi*dz/2*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=MFINTTE; 

dv=r*dr*pi*dz/2; 

C=denss*cs*dv; 

CalculateNode(T[i+l][j],T[i- 1 ] D']>0,T[i][j- 1 ]); 

}else if(j>p3){ /* steel nodes */ 

Rip=dr/((r+dr/2)*pi*dz*ks); 

Rim=dr/ ((r-dr/2) *pi*dz*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=dz/(r*dr*pi*ks); 

dv=r*dr*pi*dz; 

C=denss*cs*dv; 

CalculateNode(T[i+ 1 ]0],T[i- 1 ] [j] ,T[i] [j+ 1 ] ,T[i] [j- 1 ]); 
}else{ _ /* nodes on boundary */ 

Rip=dr/((r+dr/2)*pi*dz*ks); 
Rim=dr/((r-dr/2)*pi*dz*kt); 
Rjm=dz/((pi/2*(dr*dr/4+r*dr)*ks)+(pi/2*(r*dr- 
dr*dr/4)*kt)); 

Rjp=dz/((pi/2*(dr*dr/4+r*dr)*ks)+(pi/2 lt: (r*dr- 

dr*dr/4)*ks)); 

C=denss*cs*pi/2*(r*dr+dr*dr/4)*dz + 
denst*ct*pi/2*(r*dr-dr*dr/4)*dz; 
CalculateNode(T[i+l]D],T[i-l][j],T[i]0+l],T[i]0-l]); 

} 

}else if(j==pl){ /* nodes along outer radial boundary */ 
if(i>p2){ /* nodes along boundary */ 

Rip=dr/((r+dr/2)*pi*dz/2)*( 1 /kt+ 1 /ks); 
Rim=dr/((r-dr/2)*pi*dz/2)*( 1 /kt+ 1/ks); 

Rjm=dz/(r *dr*pi *kt) ; 

Rjp=dz/(r*dr*pi*ks); 

dv=r*dr*pi*dz; 

C=.5*denst*ct*dv + .5*denss*cs*dv; 
CalculateNode(T[i+l ] [j],T[i- 1 ] [j] ,T[i] [j+ 1 3 >T[i] [j- 1 ]); 
}else{ /* nodes in target */ 

Rip=dr/((ri-dr/2)*pi*dz*kt); 

Rim=dr/((r-dr/2)*pi*dz*kt); 
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Rjm=dz/(r*dr*pi*kt); 

Rjp=dz/(r*dr*pi*kt); 

dv=r*dr*pi*dz; 

C=denst*ct*dv; 

^ CaiculateNode(T[i+l][j],T[i-l][j],T[i]0+l],T[i][j-l]); 

}else if(j==p3){ /* nodes along inner radial boundary */ 
if(i<p2){ /* nodes along boundary */ 

Rip=dr/((r+dr/2)*pi*dz/2)*( 1 /kt+1 /ks); 
Rim=dr/((r-dr/2)*pi*dz/2)*(l/kt+l/ks); 
Rjm=dz/(r*dr*pi*kt); 

Rjp=dz/ (r * dr *pi *ks) ; 
dv=r*dr*pi*dz; 

C= 5 *denst*ct*dv + .5*denss*cs*dv; 

CalculateNode(T[i+l]D],T[i-l]D],T[i]0+l],T[i]D'-l]); 

}eise{ /* nodes in steel */ 

Rip=dr/((r+dr/2)*pi*dz*ks); 
Rim=dr/((r-dr/2)*pi*dz*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=dz/(r*dr*pi*ks); 

dv=r*dr*pi*dz; 

C=denss*cs*dv; 

^ CaicuiateNode(T[i+l][j],T[i-l][j],T[i][j+l],T[i]|j-l]); 

}else if(j==NUMZ-l){ /* nodes along insulated boundary */ 
Rip=dr/((r+-dr/2)*pi*dz/2*ks); 
Rim=dr/((r-dr/2)*pi*dz/2*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=INFINITE; 

dv=r*dr*pi*(dz/2); 

C=denss*es*dv; 

CalculateNode(T[i+l]0],T[i-l]0],0,T[i][j-l]); 

}else if(j<pl){ /* nodes in target front interior */ 
Rip=dr/((rfdr/2)*pi*dz*kt); 
Rim=dr/((r-dr/2)*pi*dz*kt); 

Rjm=dz/(r*dr*pi*kt); 

Rjp=dz/(r*dr*pi*kt); 

dv=T*dr*pi*dz; 

C=denst*ct*dv; 

CalculateNode(T[i+l][j],T[i-l]lj],T[i][j+l],T[i][j-l]); 

}else if(j<p3){ /* nodes in target stem region */ 

if(i<p2){ /* target nodes */ 

Rip=dr/((r+dr/2)*pi*dz*kt); 


Rim=dr/ ((r-dr/2) *pi * dz * kt) ; 

Rjm=dz/(r*dr*pi*kt); 

Rjp=dz/(r*dr*pi*kt); 
dv=r * dr *pi* dz; 

C=denst*ct*dv; 

CalculateNode(T[i+ 1 ][j] ,T[i- 1 ] 0] ,T[i] [j+ 1 ] ,T[i] [j- 1 ]); 
}else{ /* steel nodes */ 

Rip=dr/((r+dr/2)*pi*dz*ks); 
Rim=dr/((r-dr/2)*pi*dz*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=dz/ (r*dr*pi*ks) ; 
dv=r*dr*pi*dz; 

C=denss*cs*dv; 

CalculateNode(T[i+l][j],T[i-l]0],T[i]0+l],T[i][j-l]); 

} 

}else{ /* remaining steel nodes */ 

Rip=dr/((r+-dr/2)*pi*dz*ks); 
Rim=dr/((r-dr/2)*pi*dz*ks); 

Rjm=dz/(r*dr*pi*ks); 

Rjp=dz/(r*dr*pi*ks); 

dv=r !!: dr*pi*dz; 

C=denss*cs*dv; 

CalcuIateNode(T[i+l]0],T[i-l]D],T[i]0+l],T[i]0-l]); 

} /* end else */ 

} /* end i loop */ 

}while(!converged); /* end inner do loop */ 
gotoxy(l,5); 

printf("current time: %4.11f ms\n",t*1000); 

printf("temp at center of target is %4.01f %4.01f %4.01f %4.01f %4.01f\n", 

T[0] [0] S T[0] [ 1 J,T[0][2],T[0] [3] 5 T[0] [4]); 
printf("temp at edge of target is %4.01f %4.01f %4.0if %4.01f %4.01f\n", 
T[NUMR-l][0],T[NUMR-l][l],Tp4UMR-l][2], 
T[NUMR- 1 ] [3],T[NUMR- 1 ] [4]); 
printf("\nfront stem insulated\n"); 

for(l=NUMR- 1 ;1>=0:1 — ) 

printf("%4.01f %4.01f %4.01f %4.0!f %4.0!f %4.01f %4.01f\n" ) 

T[l] [0],T[1] [p 1 /2],T[l][p 1 ],T[1] [(p3+p 1 )/2],T[l] [p3], 

T[1][(NUMZ-1 + P 3)/2],T[1][NUMZ-1]); 


fprintf(Front, "\n%4. 1 If ”,t* 1 000); 


fprmtf(Side,"\n%4. 1 If ",t*1000); 
fprintf( C enter, "\n%4 . 1 If ",t* 1 000); 
for(l=0;l<NUMZ;l++){ 

fprintf(Side,"%4.01f ", 1 ,8*T[NUMR- 1 ][1]); 
fprintf(Center,"%4.01f ", 1 .8*T[0][1]); 

} 

for(l=0;l<NUMR;l++) fprmtf(Front,"%4.01f ",1.8*T[l][0j); 

for(l=0;l<NUMR;l++) { 

CheckF orCritical(T[l] [0] ,T critt,l,0); 
if(l<p2) CheckForCritical(T[l][p3],Tcrits,l,p3); 
else CheckForCritical(T[lJ [p l],Tcrits,l,p 1 ); 

} 

for(l=0;l<=p3 ;1++) { 

if(l<p 1 ) CheckForCritical(T[NUMR- 1 ] [l],Tcrits,NUMR- 1 ,1); 
else CheckForCritical(T[p2][l],Tcrits,p2,l); 

} 

CheckForCritical(T[NUMR- 1 ] [p 1 ], Tents, NUMR- 1 ,p 1 ); 
CheckForCritical(T[p 1 ] [p2],Tcrits,p 1 ,p2); 


for(i=0;i<NUMR;i++) 

for(j=Oii<NUMZy++)Told[i][i]=T[i]D]; 


}while(l); 

_setcursortype(_NORMALCURSOR); 

fclose(Front); 

fclose(Side); 

fclose(Center); 

} 


void CheckForCritical(double temp, double tcrit, int noder, int nodez){ 
if(temp>=terit) { 
gotoxy(l,20); 

printf("\n\n\nsafe residence time is %4.11f ms", t* 1000); 

printf(" at r=%1.2e and z=%1.2e",noder*dr*39.37,nodez*dz*39.37); 

printf("\nT=%4,01f Tcrit=%4.0Ifn", temp, tcrit); 

exit(0); 

} 

} 

void CalculateNode(double Tip, double Tim, double Tjp, double Tjm){ 


SumR=l/Rip + 1/Rim + 1/Rjp + 1/Rjm; 
SumTR=Tip/Rip + Tim/Rim + Tjp/Rjp + Tjm/Rjm; 
T[i]0]=(SumTR + (C/dt)*Told[i][j])/(SumR + C/dt); 
if(fabs(T[i][j]-To) < EPS){ 
converged*=l ; 

}else converged*=0; 
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Table Bl. Temperature Profiles (°R) Along Exposed Flat Surface 
of Aluminum Target (Figure 13) 


Time 

ms 

Centerline 

0.05 inch 

0.1 inch 

0.15 inch 

0.2 inch 

Outer 

Edge 

1 

595 

595 

595 

595 

596 

607 

10 

805 

805 

806 

808 

820 

861 

20 

918 

918 

921 

929 

951 

997 

30 

1002 

1003 

1009 

1023 

1051 

1098 

40 

1074 

1076 

1084 

1102 

1134 

1181 

51 

1144 

1148 

1158 

1179 

1213 

1260 


Table B2. Temperature Profiles (°R) Along the Exposed Flat 
Surface of Copper Target (Figure 14) 

Time 

Centerline 

0.05 inch 

0.1 inch 

0.15 inch 

0.2 inch 

Outer 

ms 






Edge 

i 

576 

576 

576 

576 

576 

584 

50 

929 

931 

939 

953 

974 

1003 

100 

1114 

1118 

1127 

1143 

1165 

1194 

150 

i a/* i 

IZ04 

1268 

1277 

1292 

1314 

1341 

200 

1393 

1396 

1405 

1420 

1440 

1466 

250 

1507 

1510 

1519 

1533 

1552 

1577 

300 

1611 

1614 

1622 

1635 

1654 

1678 

350 

1706 

1709 

1717 

1730 

1748 

1771 

384 

1768 

1770 

1778 

1791 

1809 

1831 


Table B3. Temperature Profiles (°R) Along Centerline of 
Aluminum 1 arget and ibteel fielder (Figure 15) 


Time 

ms 

0 

inch 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

1 

2 

3 

inch 

1 

595 

540 

540 

540 

540 

540 

540 

540 

540 

540 

10 

805 

550 

540 

540 

540 

540 

540 

540 

540 

540 

20 

918 

584 

543 

541 

540 

540 

540 

540 

540 

540 

30 

1002 

628 

551 

543 

542 

541 

540 

540 

540 

540 

40 

1074 

676 

566 

549 

547 

543 

540 

540 

540 

540 

51 

1144 

730 

589 

560 

554 

547 

541 

540 

540 

540 


52 

I 


i 




i 
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Table B4. Temperature Profiles (°R) Along Centerline of 
Copper Target and Steel Holder (Figure 16) 


Time 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 1 

2 

3 

ms 

inch 









inch 

i 

576 

540 

540 

540 

540 

540 

540 

540 540 

540 

540 

i 50 

929 

676 

581 

558 

553 

547 

542 

540 540 

540 

540 


1114 

842 

694 

628 

599 

576 

555 

546 542 

540 

540 

150 

1264 

988 

813 

716 

66 3 

620 

584 

565 550 

541 

541 


1393 

1116 

927 

809 

735 

676 

626 

596 565 

545 

545 

250 

1507 

1232 

1034 

902 

813 

740 

677 

636 589 

552 

554 


1611 

1338 

1135 

993 

893 

809 

735 

684 619 

565 

567 

J 350 

1706 

1437 

1231 

1082 

974 

881 

797 

737 655 

582 

585 

384 

1768 

1501 

1294 

1141 

1029 

931 

841 

776 682 

597 

600 

1 

Table B5. Temperature Profiles (°R) Along Exposed Curved Surface 

i 

i 

1 


of Aluminum Target and Steel Holder (Figure 17) 



Time 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

1 

2 

3 












) 1 

ESS 

552 

552 

552 

552 

552 

550 

550 

550 


t A 

JLU 

££ 1 
uu x 

£A0 

can 

~>yy 

rnn 

jyy 

mn 

jyy 

634 

>1 o 

Old 

/"i n 

Old 

Old 

Jr i 

016 

20 

997 


630 

628 

628 

684 

669 

669 

669 

669 

30 


733 

658 

651 

651 

717 





40 

1181 

794 

688 

672 

671 

744 

739 

739 

739 

739 

51 


858 

722 

694 

691 

768 

768 

768 

768 

768 

i 

Table B6. Temperature Profiles (°R) Along Exposed Curved Surface 




of Copper Target and Steel Holder (Figure 18) 



Time 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

1 

2 

3 

ms 

inch 









inch 

1 

584 

548 

548 

548 

548 

550 

550 

550 

550 

550 

50 

1003 

755 

662 

640 

638 

695 

765 

765 

765 

765 


1194 

927 

783 

720 

701 

760 

866 

867 

867 

867 

150 

1341 

1070 


808 

770 

821 

942 

944 

944 

944 

200 

1466 

1195 

1011 

900 

844 

887 

1006 

1009 

1009 

1009 

250 

1577 

1308 

1116 

991 

922 

956 

1064 

1065 

1066 

1066 


1678 

1411 

1214 

1080 

1001 

1029 

1119 

1117 

1118 

1118 

350 

1771 

1508 

1308 

1167 

1081 

1102 

1174 

1165 

1166 

1166 

384 

1831 

1569 

1369 

1225 

1135 

1152 

1210 

1195 

1197 

1197 


1 


p 
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/* This program calculates the velocity of particles in the rocket exhaust 

and on impact with the probe. It also calculates gas velocity, Mach number, 
and density in the exhaust flow */ 

#include <stdio.h> 

#include <math.h> 

#include <conio.h> 

#include <stdlib.h> 

#define N 1000 

FILE * density, *mach, * gas vel, * AR, * impct, * dragcoefF; 

FILE *pl0,*p50; 

FILE *pl00,*p300,*p500; 

FILE *pl000,*p3000,*p5000; 

double M,IVp,Vp_iplus 1 ,Vg_i, Vg_iplus 1 ,RHOg_i,RHOg_iplus 1 ; 

double Vpl0,Vp50; 

double Vp 1 00, Vp3 00, Vp500; 

double Vp 1 000, Vp3 000, Vp5 000; 

double RHOp,Cd,dx,pi; 

void main(){ 
double a,b,c; 

double Astar,gamma,Ac,x,xstar,xe,RHOo,To,T,R,as,Ae,Po; 

double Ar,f,dfdm; 

intij; 

void Open_Files(void); 

void Close_Files(void); 

double Particle V elocity(double, double): 

double ImpactV elocity(double, double); 

double interp(double, double, double,double, double); 

pi=4.0*atan(1.0);gamma=1.2; 

xe=48.0;xstar=xe/3.0;dx=xe/(N+l); 

Po=473.0*144;R=2093.0;To=6300.0; 

RHOo=Po/(R*To)*32. 1 74; 

Ac=4.0;Ae=6,0;Astar=l .0; 

RHOp=l 15.4879; 

clrscr(); 

Open_Files(); 
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T=To/(l .0+(gamma-l .0)/2.0); 
as=sqrt(gamma*R*T); 

for(x=0.0;x<=xe;x+=dx){ 

/* calculate Area ratio */ 
if(x/xstaK=l){ 

// Ar = ( 1. 0- Ac/Astar) *x/xstar + Ac/Astar; 

Ar=((Ac-.5*(Ac+Astar))*cos(pi*x/xstar)+.5*(Ac+Astar))/Astar; 

M = 0.0; 

}else{ 

// Ar = (Ae/Astar - 1.0)*(x-xstar)/(xe-xstar) + 1.0; 

Ar=((Ae-.5*(Ae+Astar))*cos(pi*(x-xe)/(xstar-xe))+.5*(Ae+Astar))/Astar; 

M = 1.5; 

} 

/* calculate M */ 

a=pow((gamma+ 1 ,0)/2.0,(gamma+ 1 .0)/(2.0*(gamma- 1 .0))); 
b=(gamma- 1 .0)/2 .0; 
c=-(gamma+l .0)/(2.0*(gamma- 1 .0)); 
do{ 

f=a*M*pow((l .0+b*M*M),c)-l .G/Ar; 

' dfdm = a*M*c*pbw((l .0+b*M*M),c-l .0)*2.0*b*M + a*pow((l .0+b*M*M),c); 

ivx—ivx-iy uxuiii, 

} while(fabs(f/dfdm)>.0000 1 ); 
if(x=0){ 

Vg_i = as*sqrt( (gamma+1.0)/2.0 * M*M / (L0+(gamma-1.0)/2.0*M*M) 

RHOgJ = RHOo*pow(( 1 .0+(gamma- 1 .0)/2.0*M*M),-l .G/(gamina- 1 .0)); 

Vp 1 0=Vp50=0;//V g_i; 

Vp 1 00=Vp3 00=Vp500=0;// VgJ; 

Vp 1 000=Vp3000=Vp5000=0;//Vg_i; 

}else{ 

Vg_iplusl = as*sqrt( (gamma+1.0)/2.0 * M*M / (L0+(gamma-l ,0)/2.0*M*M) ); 
RHOg_iplusl = RHOo*pow((l .0+(gamma- 1 .0)/2.0*M*M),- 1 .0/(gamma- 1 .0)); 

if(M<0.2)Cd=0.48; 

else if(M<=0.4)Cd=interp(0.2,0.48,0.4,0.5,M); 

else if(M<=0.6)Cd=interp(0.4,0.5, 0.6, 0.55, M); 
else if(M<=0.8)Cd=interp(0.6, 0.55, 0.8, 0.62, M); 
else if(M<=1.0)Cd=interp(0.8, 0.62, 1.0, 0.8, M); 
else if(M<=1.2)Cd=interp(l. 0,0.8, 1 .2,0.93 ,M); 
else if(M<=l .4)Cd=interp(l. 2, 0.93,1 .4, 1 .05, M); 


else if(M<=2.0)Cd=interp(l .4,1 .05,2.0, 1 .04, M); 

else if(M<=3 .0)Cd=interp(2.0, 1 .04,3. 0,0.95, M); 
else Cd=.95; 

Vpjiplus 1 =ParticleVelocity(0.0000033 * 1 0,Vp 1 0); 
fprintf(pl0,"%5.31f %5.01f\n",x/xe,Vpl0); 

VplO=Vp_iplusl; 

Vp_iplusl=ParticleVelocity(0.0000033*50,Vp50); 
fprintf(p50,"%5.31f %5.01f\n",x/xe,Vp50); 

Vp50=Vp_iplus 1 ; 

Vp Jplus 1 =ParticleVelocity(0. 0000033* 1 00, Vp 1 00); 
fprintf(pl00,"%5.31f %5.01f\n",x/xe,Vpl00); 

Vp 1 00=Vp_iplus 1 ; 

Vp_iplusl=ParticleVelocity(0.0000033*300,Vp300); 

fprintf(p300,"%5.31f %5.01f\n",x/xe,Vp300); 

Vp3 00=Vp_iplus 1 ; 

Vp_iplusl=ParticleVelocity(0.0000033*500,Vp500); 
fprintf(p500,"%5.31f %5.01f\n",x/xe,Vp500); 

V p5 00= V p_iplus 1 ; 

Vp_iplusl=ParticleVelocity(0.0000033*1000,Vpl000); 

fprintf(pl000, ,,o /o5.31f %5.01f\n",x/xe,Vpl000); 

Vp 1 000=Vp_iplus 1 ; 

Vp_iplusl=ParticleVelocity(0.0000033*3000,Vp3000); 
IVp=ImpactVelocity(0.0000033*3000,Vp3000); 
fpriiitf(p3000,"%5.31f %5.01f\n",x/xe,Vp30G0); 

Vp3 000=Vp_iplus 1 ; 

Vp_iplusl=PartieleVeloeity(0.0000033*5000,Vp5000); 
fprintf(p5000,"%5.31f %5.01f\n",x/xe,Vp5000); 
Vp5000=Vp_iplusl ; 

Vg_i=Vg_iplusl; 

RHOg_i=RHOgiplus 1 ; 

} 


fprintf(density, "%5 .31f %lf\n",x/xe,RHOg_i); 
fbrintiTAR. "%5 ,3If %5.31f\n",x/xe,Ar); 


fprmtf(mach,"%5.31f %5.31f\n",x/xe,M); 
fprintf(gasvel," 0 /o5.31f %5.01f\n”,x/xe,Vg_i); 
fprintf(dragcoeff,"%5.31f %5.31f\n",M,Cd); 

} 

x-=dx; 

IVp=ImpactVelocity(0. 0000033*. 00 1 ,Vg_i); 
fprintf(impct, "%5 ,31f %5.01f %5.01f\n",.00 1 ,Vg_i,IVp); 

rVp=ImpactVelocity(0.0000033*.0 1 ,Vg_i); 
fprmtf(impct,"%5.31f %5.01f %5.01f\n",.01,Vg_i,IVp); 

rVp=ImpactV elocity(0.0000033 * . 1 , Vg_i); 
fprintf(impct,"%5.31f %5.01f %5.01f\n",.l,Vg_i,IVp); 

IVp=ImpactV elocity(0. 000003 3*1 0, Vp 1 0); 
fprintf(impct,"%5.31f %5.01f %5.01f\n",10.0,Vpl0IVp); 

IVp=ImpactVelocity(0.0000033*50,Vp50); 
fprmtf(irnpct,"%5.3If %5.01f %5.01f\n",50.0,Vp50,IVp); 

IVp=ImpactVelocity(0.0000033 * 1 00, Vp 1 00); 
fprintf(impct, "%5 .3 If %5.01f %5.01f\n",100.0,Vpl00,IVp); 

rVp=ImpactVelocity(0.0000033*300.Vp300); 
fprintf(impct, "%5 .3 If %5.01f %5.01f\n", 300.0, Vp300,IVp); 

IVp=ImpactVelocity(0.0000033*500,Vp500); 
^rintf(impct,"%5.31f %5.01f %5.01f\n", 500.0, Vp500, IVp); 

IVp=ImpactVelocity(0.0000033*l 000, Vp 1 000); 
fprin1f(impct,"%5.31f %5.01f %5.01f\n", 1 000.0, Vpl 000, IVp); 

IVp=ImpactVelocity(0.0000033*3000,Vp3000); 
fprintf(impci ,,0 /o5.31f %5.01f %5.01f\n",3000.0,Vp3000,IVp); 

IVp=ImpactVelocity(0.0000033*5000,Vp5000); 
fprmtiTpSOOO, "°/o5 .3 If %5.0If %5.01f\n", 5000.0, Vp5000, IVp); 

Close_Files(); 

printf( "Calculations completed successfully^"); 

} 


void Open_FiIes(){ 


if ((density = fopen("c:\\partvel\\density.dat", "wt")) — NULL){ 

fprintf(stderr, "Cannot open output file DENSITY.DAT.\n"); 
exit(l); 

} 

if ((mach = fopen("c:\\partvel\\mach.dat" } "wt")) = NULL){ 

lprintf(stderr, "Cannot open output file MACH.DAT.\n"); 
exit(l); 

} 

if ((gasvel = fopen("c:\\partvel\\gasvel.dat", "wt")) = NULL){ 

lprintf(stderr, "Cannot open output file GASVEL.DAT.\n"); 
exit(l); 

} 

if ((AR = fopen("c:\\partveI\\AR.dat", "wt")) = NULL){ 

fprintf(stderr, "Cannot open output file AR.DAT.\n"); 
exit(l); 

} 

if ((plO = fopen("c:\\partvel\\plOe.dat" s "wt")) = NULL){ 

fprintf(stderr, "Cannot open output file P10.DAT.\n"); 
exit(l); 

} 

if ((p50 = fopen("c:\\partvel\\p50e.dat", "wt")) = NULL){ 

fprintf(stderr, "Cannot open output file P50.DAT.\n"); 
exit(l); 

} 

if ((pi 00 = fopen("c:\\partveI\\p 100e.dat", "wt")) == NULL){ 
fprintf(stderr, "Cannot open output file P100.DAT.\n"); 
exit(l); 

} 

if ((p300 = fopen("c:\\partvel\\p300e.dat", "wt")) = NULL){ 
fprintf(stderr, "Cannot open output file P300.DAT.\n"); 
exit(l); 

} 

if((p500 = fopen("c:\\partvel\\p500e.dat", "wt")) = NULL){ 
fprintf(stderr, "Cannot open output file P500.DAT.\n"); 
exit(l); 

} 

if ((plOOO = fopen("c:\\partvel\\p 1 000e.dat", "wt")) = NUJLL){ 
ft)rintf(stderr, "Cannot open output file PiOQO.DATAn"): 
exit(l); 

} 

if ((p3000 = fopen("c:\\partvel\\p3000e.dat", "wt")) = NULL){ 
fi3rintf(stderr, "Cannot open output file P3000.DAT.\n"); 
exit(l); 


} 

if (CpSOOO = fopen("c:\\partvel\\p5 OOOe.dat", "wt")) = NULL){ 
fprintf(stderr, "Carnot open output file P5000.DAT.\n"); 
exit(l); 

} 

if ((impct = fopen("c:\\partveI\\impvele.dat", "wt")) = NULL){ 
fprmtf(stderr, "Cannot open output file P5000.DAT.\n"); 
exit(l); 

} 

if ((dragcoeff = fopen("c:\\partvel\\cd.dat", "wt")) — NU LL ){ 
fprintf(stderr, "Cannot open output file cd.DAT.\n"); 
exit(l); 

} 

} 

void Close_Files(){ 
fclose(density); 
fclose(AR); 
fclose(mach); 
fclose(gasvel); 
fclose(plO); 
fclose(p50); 

iwius>e*jpiuuj 3 

fclose(p300); 

fclose(p500); 

fclose(pl000); 

fclose(p3000); 

fclose(p5000); 

fclose(impct); 

fclose(dragcoeff); 

} 

double ParticleV elocity(double D, double Vpi){ 
double a,b,c,d,A,B,C,srt; 
double Vgbar,RHObar; 

V gbar=.5*(V: g_iplus 1 + V g_i); 

RHObar=.5 *(RHOg_i+RHOg_ipIus 1 ); 

a=2.0*RHOp*D/(Cd*RHObar); 

b=-1.0; 

c=2.0*Vgbar; 

d=- V gbar* V gbar; 


A=a/dx + b/2.0; 

B=b*Vpi + c; 

C=Vpi*Vpi*(b/2.0-a/dx) + c*Vpi + 2.0*d; 

srt=sqrt(B *B-4 .0 * A * C); 
if((-B+srt)<0)srt=-srt; 

return (-B+srt)/(2.0*A); 

} 

double ImpactVelocity(double D, double Vp){ 
double del; 

del=.143*exp(3.24/(M*M))*.25/12; 

return Vp*exp(-Cd*del*RHOg_i/(4.0/3 .0*D*RHOp)); 

} 

double interp(double xl, double fl, double x2, double £2, double x){ 
return fl +(£2-fl )/(x2-x 1 )*(x-x 1 ); 

} 
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Table Dl. Nozzle Flow Properties as a Function of X/X e 


X/Xe 

A/A* 

Mach Number 

Gas Velocity 
ft/sec 

Gas Density 
lbm/ft 3 

Chamber 

4 

0.15 

595 

0.1643 

0.1 

3.383 

0.178 

707 

0.1636 

0.2 

2.039 

0.306 

1210 

0.1587 

0.3 

1.075 

0.737 

2854 

0.1037 

Throat 

1 

1 

3773 

0.1276 

0.4 

1.121 

1.383 

5040 

0.0693 

0.5 

1.736 

1.922 

6534 

0.0345 

0.6 

2.732 

2.322 

7445 

0.0192 

0.7 

3.895 

2.599 

7987 

0.0126 

0.8 

4.971 

2.781 

8307 

0.0095 

0.9 

5.728 

2.884 

8476 

0.0081 

Exit 

6 

2.917 

8529 

0.0076 


Table D2. Velocity (ft/sec) of a 10 pm Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/X e 

X e -1 inch 

X e =6 inch 

X e =12 inch 

X e =24 inch 

X e =48 inch 

Chamber 

0 

0 

0 

0 

0 

0.1 

439 

579 

612 

636 

654 

0.2 

684 

907 

974 

1029 

1073 

0.3 

1460 

2039 

2220 

3172 

2488 

Throat 

1961 

2740 

2979 

2369 

3324 

0.4 

2977 

3975 

4245 

4453 

4609 

0.5 

3998 

5210 

5538 

5795 

5990 

0.6 

4587 

5949 

6321 

6613 

6835 

0.7 

4956 

6413 

6811 

7122 

7359 

0.8 

5214 

6730 

7141 

7459 

7699 

0.9 

5410 

6961 

7376 

7692 

7927 

Exit 

5573 

7138 

7548 

7855 

8076 
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Table D3. Velocity (ft/sec) of a 50 Jim Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/Xe 

X e =l inch 

X e =6 inch 

X e =12 inch 

X e =24 inch 

X e =48 inch 

Chamber 

0 

0 

0 

0 

0 

0.1 

279 

456 

517 

566 

603 

0.2 

451 

710 

801 

883 

954 

0.3 

926 

1523 

1758 

1975 

2165 

Throat 

1235 

2046 

2363 

2654 

2907 

0.4 

1925 

3094 

3513 

3873 

4165 

0.5 

2653 

4143 

4654 

5088 

5441 

0.6 

3068 

4749 

5322 

5811 

6210 

0.7 

3328 

5130 

5743 

6265 

6693 

0.8 

3510 

5395 

6033 

6577 

7019 

0.9 

3652 

5597 

6251 

6806 

7253 

Exit 

3774 

5762 

6426 

6983 

7428 


Table D 4. Velocity (ft/sec) of a 100 Particle in the Exhaust Flow 

of Various Size Rockets as a Function of X/X e 


X/Xe 

X e =l inch 

X e =6 inch 

Xe=12 inch 

X e =24 inch 

X e =48 inch 

Chamber 

0 

0 

0 

0 

0 

0.1 

218 

388 

456 

517 

566 

0.2 

359 

611 

710 

801 

883 

0.3 

728 

1284 

1523 

1758 

1975 

Throat 

968 

1722 

2046 

2363 

2654 

0.4 

1520 

2641 

3094 

3513 

3873 

0.5 

2112 

3577 

4143 

4654 

5088 

0.6 

2451 

4113 

4749 

5322 

5811 

0.7 

2662 

4450 

5130 

5743 

6265 

0.8 

2812 

4685 

5395 

6033 

6577 

0.9 

2928 

4865 

5597 

6251 

6806 

Exit 

3029 

5017 

5762 

6426 

6983 
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Table D5. Velocity (ft/sec) of a 300 (j.m Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/X e 

X e =l inch 

Xe=6 inch 

Xe=12 inch 

X e =24 inch 

X e =48 inch 

Chamber 


0 

0 



0.1 


279 

347 

417 

482 

0.2 

236 

451 

551 

652 

749 

0.3 

475 

926 

1147 

1383 

1622 

Throat 

630 

1235 

1535 

1856 

2180 

0.4 

994 

1925 

2371 

2831 

3274 

0.5 

1397 

2653 

3233 

3817 

4363 

0.6 

1627 

3068 

3726 

4383 

4996 

0.7 

1772 

3328 

4035 

4738 

5394 

0.8 

1874 

3510 

4251 

4986 

5670 

0.9 

1954 

3652 

4418 

5176 

5879 

Exit 

2024 

3774 

4559 

5334 

6049 


1 able SJb. V elocity (ft/sec) of a 500 fun Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/Xe 

X e =l inch 

X e =6 inch 

X e — 12 inch 

X e =24 inch 

X e =48 inch 

Chamber 


0 

0 

0 


0.1 


233 

296 

365 

435 

0.2 


382 

477 

578 

678 

0.3 

383 

777 

982 

1208 

1446 

Throat 

507 

1035 

1311 

1618 

1942 

0.4 

803 

1621 

2040 

2492 

2951 

0.5 

1132 

2249 

2803 

3388 

3965 

0.6 

1321 

2607 

3239 

3901 

4550 

0.7 

1439 

2831 

3511 

4222 

4917 

0.8 

1523 

2989 

3703 

4447 

5173 

0.9 

1589 

3112 

3851 

4620 

5368 

Exit 

1646 

3218 

3978 

4766 

5530 
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Table D7. Velocity (ft/sec) of a 1000 urn Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/Xe 

X e =l inch 

X e =6 inch 

X e =12 inch 

Xe=24 inch 

X e =48 inch 

Chamber 

0 

0 

0 

0 

0 

0.1 

82 

179 

233 

296 

365 

0.2 

141 

297 

382 

477 

578 

0.3 

283 

601 

777 

982 

1208 

Throat 

374 

798 

1035 

1311 

1618 

0.4 

593 

1256 

1621 

2040 

2492 

0.5 

839 

1755 

2249 

2803 

3388 

0.6 

981 

2041 

2607 

3239 

3901 

0.7 

1070 

2219 

2831 

3511 

4222 

0.8 

1133 

2346 

2989 

3703 

4447 

0.9 

1182 

2445 

3112 

3851 

4620 

Exit 

1226 

2530 

3218 

3978 

4766 


Table DS. Veiocity (ft/sec) of a 3000 Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/Xe 

X e =l inch 

X e =6 inch 

Xe=12 inch 

X e =24 inch 

Xe=48 inch 

Chamber 

0 

0 




0.1 

49 

112 

151 

200 

258 

0.2 

86 

191 

254 

331 

421 

0.3 

171 

383 

512 

670 

859 

Throat 

226 

507 

678 

891 

1145 

0.4 

359 

803 

1071 

1400 

1789 

0.5 

510 

1132 

1502 

1951 

2473 

0.6 

597 

1321 

1748 

2266 

2863 

0.7 

651 

1439 

1903 

2463 

3107 

0.8 

690 

1523 

2012 

2602 

3279 

0.9 

721 

1589 

2098 

2711 

3413 

Exit 

748 

1646 

2173 

2805 

3527 
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Table D9. Velocity (ft/sec) of a 5000 pm Particle in the Exhaust Flow 
of Various Size Rockets as a Function of X/X e 


X/X e 

X e =l inch 

X e =6 inch 

X e =12 inch 

X e =24 inch 

X e =48 inch 

Chamber 

0 

0 

0 

0 

0 

0.1 

39 

90 

122 

163 

214 

0.2 

67 

153 

206 

273 

354 

0.3 

135 

307 

414 

550 

717 

Throat 

177 

406 

549 

730 

954 

0.4 

282 

643 

868 

1151 

1497 

0.5 

401 

909 

1221 

1612 

2082 

, 0.6 

470 

1062 

1424 

1876 

2416 

0.7 

513 

1158 

1552 

2041 

2625 

0.8 

544 

1226 

1642 

2158 

2773 

0.9 

568 

1280 

1713 

2250 

2888 

Exit 

589 

1327 

1774 

2329 

2987 


Table DIO. Particle Velocity at the Nozzle Exit Plane, V pe (ft/sec), 
and Particle Impact Velocity, V p (ft/sec) 


Dp 

fim 

X e =l 

Vpe 

inch 

v P 

X e =6 

Vpe 

inch 

v P 

X e =12 inch 
Vpe V p 

X e =24 inch 
Vpe V p 

X e =48 inch 
Vp e V p 

0.001 

8529 

0 

8529 

0 

8529 

0 

8529 

0 

8529 

0 

0.01 

8529 

16 

8529 

16 

8529 

16 

8529 

16 

8529 

16 

0.1 

8529 

4551 

8529 

4551 

8529 

4551 

8529 

4551 

8529 

4551 

10 

5574 

5539 

7140 

7095 

7549 

7502 

7856 

7807 

8078 

8027 

50 

3775 

3770 

! 5764 

5757 

6427 

6419 

6985 

6976 

7429 

7420 

100 

3030 

3028 

5018 

5015 

5764 

5760 

6427 

6423 

6985 

6980 

300 

2024 

2024 

3775 

3774 

4560 

4559 

5335 

5334 

6051 

6050 

500 

1647 

1647 

3219 

3219 

3979 

3979 

4767 

4767 

5531 

5530 

1000 

1226 

1226 

2531 

2531 

3219 

3219 

3979 

3979 

4767 

4767 

3000 

748 

748 

1647 

1647 

2173 

2173 

2806 

2806 

3528 

3528 
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